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All the leading-twist parton distribution functions are calculated in a spectator model of the 
nucleon, using scalar and axial-vector diquarks. Single gluon rescattering is used to generate T-odd 
distribution functions. Different choices for the diquark polarization states are considered, as well 
as a few options for the form factor at the nucleon-quark-diquark vertex. The results are listed 
in analytic form and interpreted in terms of light-cone wave functions. The model parameters 
are fixed by reproducing the phenomenological parametrization of unpolarized and helicity parton 
. . . distributions at the lowest available scale. Predictions for the other parton densities are given and, 

OO ' whenever possible, compared with available phenomenological parametrizations. 

o ■ 

O . PACS numbers: 12.39.-x, 13.60.-r, IS.i 



I. INTRODUCTION 

Partonic transverse-momentum distributions (TMDs) — also called unintegrated PDFs — describe the probability 
<^ ■ to find in a hadron a parton with longitudinal momentum fraction x and transverse momentum p^. with respect to 
C^' the direction of the parent hadron momentum [1!]. They give a three-dimensional view of the parton distribution in 

momentum space, complementary to what can be obtained through generalized parton distributions [3, |3|, y, |5|, |6| . 
(j_) ' In the last years a lot of theoretical and experimental activity related to TMDs has takenplace. Crucial steps were 

tS^ , made in the understanding of factorization theorems involving TMDs {k^. factorization) [7, 8]. Some of the properties 
of TMDs have been investigated from the theoretical standpoint. For instance, positivity bounds were presented in 
Ref. 0] . Relations among these functions in the large- A'c limit of QCD were put forward in Ref. [l3| • Their behavior 
J> ■ at large x was studied in Ref. [Il[, and at high transverse momentum in Ref. [I4I. Last but not least, it was also 
C^ ' demonstrated [l3i that TMDs that are odd under naive time- reversal transformations (for brevity, T-odd) can be 
CN nonzero and must be included in the complete list of leading-twist TMDs (see, e.g., Ref. [ij, llal)- Their universality 

^^ properties are different from the standard PDFs [lg| . 

^-% : In the meanwhile, several azimuthal asymmetries were measured in semi-inclusive deep inelastic scattering (SIDIS) 

t ' and elsewhere (see Ref. 11?'] and references therein), and more experimental measurements are planned. However, not 

^D ' much phenomenological information concerning TMDs is available as yet (see, e.g., Ref. jlgj and references therein). 

The analysis of azimuthal spin asymmetries both in hadron-hadron collisions and in SIDIS led to the extraction of 

the Sivers function [19!], denoted as f^, a T-odd TMD that describes how the parton distribution is distorted by 

the transverse polarization of the parent hadron (see Ref. [20| for a comparison of various parametrizations). A 

recent attempt to extract the T-odd Boer- Mulders function, h^ [2l|, a T-odd TMD describing the distribution of 

, ; , transversely polarized partons in an unpolarized hadron, was presented in Ref. J22| . All of the above studies assume 

C^ ' a flavor-independent Gaussian distribution of the transverse momentum, although there is no compelling reason for 

this choice. 

In this context, building a relatively simple model to compute TMDs and to allow for numerical estimates is of 
great importance. From the theoretical side, this can help understanding some of the essential features of TMDs, for 
instance their relation to the orbital angular momentum of partons (see, e.g., Refs. [IJ, [23|, [2J, |2^, [2^ 123, [28|, |29J). 
From the experimental sidCja model could be useful to estimate the size of observables in different processes and 
kinematical regimes [H iH, IH, [H, Hi, HE IHl and to set up Monte Carlo simulations [13, [H, [UliO, [il . 

Although many model calculations of integrated PDFs are available, there are not so many for TMDs. In Ref. [43| 
all the leading-twist T-even functions were calculated in a spectator model with scalar and axial-vector diquarks. 
Recently, an analogous calculation has been performed in a, lig ht-cone quark model [43|. T-odd functions were 
calculated in the spectator model with scalar diquarks [ij, [30, [4J, [3 , with scalar and vector diquarks [Sy, [4y| , in the 
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MIT bag model [43, |4g| , in a constituent quark model [49'] and in the spectator model for the pion [50] . A complete 
calculation of all the leading-twist TMDs in a spectator model with scalar diquarks was presented in Ref. [23] • 

In this work, we choose a more phenomenological approach. We consider also axial- vector diquarks (in the following 
often called simply vector diquarks), necessary for a realistic flavor analysis, and we further distinguish between 
isoscalar (ud-like) and isovcctor (uu-like) spectators. We generate the relative phase necessary to produce T-odd 
structures by approximating the gauge link operator with a one gluon-exchange interaction. We consider several 
choices of form factors at the nucleon-quark-diquark vertex and several choices for the polarization states of the 
diquark. All results are presented in analytic form and interpreted also in terms of overlaps of light-cone wave 
functions, leading to a detailed analysis of the quantum numbers of the quark-diquark system. The free parameters of 
the model are fixed by reproducing the phenomenological parametrization of unpolarized and longitudinally polarized 
parton distributions at the lowest available scale. 

The paper is organized as follows. In Sec. [TTl the analytic form for all the leading-twist TMDs is discussed for 
the dipolar nucleon-diquark-quark form factor and for the light-cone choice of the diquark propagator, postponing 
the results for the other explored combinations to the Appendices |A| (T-even TMDs) and |B] (T-odd TMDs). In 
Sec. IIIH numerical results are shown and compared with phenomenological parametrizations, whenever available in 
the literature. In Sec.lIVi some conclusions are drawn. 



II. ANALYTICAL RESULTS FOR TRANSVERSE-MOMENTUM-DEPENDENT PARTON DENSITIES 

In this section we present the fundamentals of the model and we give in analytical form the results for the light-cone 
wave functions (LCWFs) and the TMDs obtained in the model. 

A. General framework 

In the following we will make use of light-cone coordinates. We introduce the light-like vectors n± satisfying 
nj_ = 0, n+ • n_ = 1, and we describe a generic 4-vector a as 

a — [a^ , a+, Uj,] (1) 

where a = a ■ rizfi. We will make use of the transverse tensor e^ = e^'^'-'Vi+^ri-j^, whose only nonzero components are 
e}^ = ~£t^ — 1- We choose a frame where the hadron momentum P has no transverse components, i.e.. 



P = 



M 
2P 



2 



- P+ 



The quark momentum can be written as 



P = 



P^ +pI 

2xP+ 



xP^,Pt 



(2) 



(3) 



In a hadronic state ]P, S) with momentum P and spin 5, the density of quarks can be defined starting from the 
quark-quark correlator (see, e.g., Ref. [l5|) 

4>(x,p,;5) = y'^|^e'^-«(P,^]V;(0)^o,a'A(OI^,^>L^_, (4) 



?+=o 



where 



U[o^^] ^ V e-'s fo dn^-M^) (5) 



is the so-called gauge link operator, or Wilson line, connecting the two different space-time points 0,^, by all possible 
ordered paths followed by the gluon field A, which couples to the quark field tp through the coupling g. The gauge 
link ensures that the matrix element of Eq. ((4|) is color-gauge invariant and arises from the interaction of the outgoing 
quark field with the spectators inside the hadron. The leading contributions of the path [0, ^] in space-time are 
selected by the hard process in which the parton distributions appear, thus breaking standard universality of the 
parton densities. For instance, in SIDIS the gauge link path in light-cone coordinates runs along 

[0,e] = (0,0,0,) -^ (0,(X3,O,) -> (0,00,00,) ^ (0,00,1,) -> (O^C^t) , (6) 



while in the Drell-Yan case it runs in the opposite direction through — oo. This fact leads to a sign difference in T-odd 
parton densities, as mentioned for the first time in Ref. [1.61]. 

Similarly to Ref. [44I, we evaluate the correlator of Eq. (j4|) in the spectator approximation, i.e. we insert a com- 
pleteness relation and at tree-level we truncate the sum over final states to a single on-shell spectator state with mass 
Mx, thus getting the analytic form 



Hx,Pt,S) 



1 



(27r)3 2(1 - x)P+ 



M^°\S) M'^"Hs) 



p^=T(X,p 



(7) 



where p is the momentum of the active quark, m its mass, and the on-shell condition (P — p)^ — M'^ for the spectator 
implies for the quark the off-shell condition 



2 / N pI + Lx(m'^) 2 

j^ = t{x,Pt) = , ^ '-+m^ 



with AI the hadron mass. 



1 -X 



L\{m'^) = xMJ^ + (1 - x)'m^ - x{\ - x)!^"^ 



(8) 
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FIG. 1: Tree-level cut diagram for the calculation of T-even leading-twist parton densities. The dashed line indicates both 
scalar and axial- vector diquarks. 

We assume the spectator to be point-like, with the quantum numbers of a diquark. Hence, the proton can couple to 
a quark and to a spectator diquark with spin (scalar X = s) or spin 1 (axial- vector X = a), as well as with isospin 
(isoscalar ud-like system) or isospin 1 (isovector uu-\ike system). Therefore, the tree- level "scattering amplitude" 
A^(o) is given by (see Fig.[T|) 



X(°)(5) = (P-p|^(0)|P,5) = 



■ysUiP,s) 



scalar diquark. 



■ e* JP ~ p, Xa) X' U{P, S) axial- vector diquark. 



(9) 



and is actually a Dirac spinor because of the understood spinorial indices of the quark field ip. The e^{P — p, Xa) is 
the 4- vector polarization of the spin-1 vector diquark with momentum P ~ p and helicity states Xa- When summing 
over all polarizations states, several choices have been used for d^'^ = J2\ ^*(\ ) ^\ 
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{P~pYn''_ + {P -pYnt 
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Ref. 
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(see 


Ref. 


[46]) 



(10) 
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The different forms for the diquark propagator correspond to different physical theories and lead to different results 
for the parton distribution functions. We have analyzed all of them except for the third one, which was extensively 
studied already in Ref. [42] ■ However, we think that the first one is preferable to the others. The motivation is 
that in the spectator model we have to take into account that the diquarks have an electric charge and can couple 
to the virtual photon in DIS. In other words, in this model the quarks are not the only charged partons in the 
proton: the diquarks are also charged partons and they have spin different from i. The scalar diquark couples only 



to longitudinally polarized photons and gives contribution to the structure function F^. This leads to a violation 
of the Callan-Gross relation, but leaves unchanged the (leading-order) interpretation of the structure function Ft 
as a charge- weighted sum of quark distribution functions. This seems the best way to reduce the phenomenological 
impact of the problem represented by the presence of the diquarks. For the vector diquark, we checked that the same 
situation occurs when only the light-cone transverse polarization states of the diquark are propagated, i.e., when the 
first choice of Eq. ()10p for the polarization sum is used. In the other cases, the diquark would give a contribution 
also to the structure function Ft- On top of this, we remark that the last choice of Eq. (I10|) for the polarization sum 
introduces unphysical polarization states of the vector diquark (see discussion in next section) . In conclusion, in the 
following we shall consider only light-cone transverse polarizations of the diquark, make only a few comments on the 
other choices, and leave the complete list of results in the Appendices. 
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FIG. 2: Tree-level cut diagram for the calculation of T-even leading- twist parton densities for an active scalar or vector diquark 
(dashed line), with a spectator quark (solid line). 

Equation ^ can be further elaborated by choosing the nucleon-quark-diquark vertex y. We choose the scalar and 
vector vertices to be 

ys^i9s{p')l, 3^^=z^^7^75, (11) 



where gx[p^) is a suitable form factor. Other choices are possible (see, e.g., Refs [3a, S^]), but we limit ourselves to 
these ones, which are the simplest. For the form factor, we explored three possible choices: 



9x{p^) = < 



g^ ■ point-like, 

2 2 

d'x' —7, TTTI dipolar, (12) 



^ ^^^Pg(p^-m^)/Ai exponential, 

where gx and Kx are appropriate coupling constants and cutoffs, respectively, to be considered as free parameters 
of the model together with the mass of the diquark Mx ■ All these parameters can in principle be different for each 
type of diquark. Only the point-like coupling can be derived from a specific Lagrangian with protons, quarks and 
diquarks as fundamental degrees of freedom, and meant to effectively describe QCD in the nonperturbative regime. 
Since our interest here is mainly phenomenological, we prefer to introduce form factors. They smoothly suppress the 
influence of high p^ — where our theory cannot be trusted — and eliminate the logarithmic divergences arising after 
Pt integration when using a point-like coupling. For later use, we note that the dipolar form factor can be usefully 
rewritten, using Eq. ([5]), as 

/ 2^ dip P^-rri^ d,p (p2 - m^) (1 - a;)2 

b^-Ajfl (p2 +ij^(Aj^)) 

In summary, we have analyzed in total nine combinations of nucleon-quark-diquark form factors and forms for the 
diquark propagator. As mentioned above, we will discuss analytical and numerical results involving the dipolar form 
factor and the first choice of Eq. pH]) (transverse diquark polarizations only), listing the formulae for the other cases 
in the Appendices [X] and [B) To keep the notation lighter, we will denote the coupling g-,^^ simply as gx from now on. 



B. Light-cone wave functions 

A convenient way to compute parton distribution functions is by making use of light-cone wave functions (LCWFs), 
as done for instance in Ref. 51]. For the scalar diquark, LCWFs can be defined as 



<"(x,p.) 



P+ "(P, \) 



(P — p)+ p^ — m? 



ys u{p, \n) , 



(14) 



where the indices \m, A^, refer to the helicity of the nuclcon and of the quark, respectively, and are constrained by 
angular momentum conservation to the "spin sum rule" \n = ^q + Lz, where Lz is the projection of the relative 
orbital angular momentum between the quark and the diquark. We use the conventions of Ref. [S^l (see also Ref. [27|). 
In standard representation, the spinors can be written as 



u{p, +) = 



n/W 



2r)+ 



u{p, 



V237v 



/ V2p+ + ni \ 

Px + iPy 

v2p^ — m 

\ Px+ iPy / 

and similarly for the nucleon spinors (changing p, tti, to P,M, respectively). We obtain 

'4>X{x,PT) = {m + xM)(l)/x {Lz = 0), 

4>^{x,Pt) = ~{px+iPy)(f>/x {Lz = +l) 

ip+ix,pT) = -[ipl{x,pT)]* {Lz = -l) 

tpZ{x,pT)^i^ti^^PT) iLz = 0), 



( -Px + iPy \ 

V2p+ + ni 

Px - iPy 

\ -V2p+ +m J 



(x,p^) = 



x{l — x) 



yr^^p2+i2(^2) 



which correspond to Eqs. (44,46) of Ref. 51]. 

For the vector diquark, LCWFs can be defined as 



(15) 



(16) 
(17) 
(18) 
(19) 

(20) 



p^ 



<\(^:Pt) - y (p _ p)+ p2 _ ^2 -M 



^^P^^sUP-pAa)yi:uiP,XN) 



(21) 



where the index Xa refers to the helicity of the vector diquark and is constrained by Ajv = Xg + Xa + Lz- The light-cone 
transverse polarization vectors are given by [5l| 



£(P-p,+) = 

e{P-p,-) = 



{P-p)^ + i{P-p) 



V2{P-p)+ 

Px - iPy 



y - '- 

' ' V2' V2 



Px + iPy 



1 i 

,0, ^, j= 



V2(l-a;)P+' ' \/2' \/2 



1 i 
0,-=, j= 



\/2(l-x)P+' ' V2' V2 



(22) 
(23) 



They satisfy the usual properties^ e(±) • e*(±) = —1, e(±) • e*(=F) — 0, and [P — p) ■ e(±) — 0. They are consistent 
with the polarization sum being expressed by the first option in Eq. (jTU]) . The LCWFs become 



ipX+ix.pT.) = — ^-(j)/x 

\ — X 

^X_{x,Pt) = ~x— ^(/./a; 

1 — a; 

V'l+(a;,pT) = (to + xM) (p/x 
V'J:_(x,pr) = 
'^++{x,Pt) = 



(i. =-1), 

(i. =+1), 

[Lz = 0), 
{Lz = +2), 
{Lz=-2), 



(24) 

(25) 

(26) 

(27) 
(28) 



^ Note that (P — p) ■ e(±) ^ 0, since e(it) do not describe transverse polarization with respect to the diquark momentum. 



^-_{x,pr) = -ijt+{^,PT) (i. =0), (29) 

V^Z+(x,p,) = ['^4_{x,p^)]* (L, = -1), (30) 

^ZA^^Pt) = [^X+{x,Pr)]* {L, = +1), (31) 

cj,{x,pl) = -^M^= ,[,,,\, , 32 

and are analogous to Eqs. (21,24) in Ref. [51|, the differences being due to the fact that here the diquark is an 
axial-vector particle rather than a vector one. Note that in our model we can only have wavefunctions with at most 
one unit of orbital angular momentum [p wave). The LCWFs with two units of orbital angular momentum (d-wave), 
V'-- and ■0++J vanish. 

If we add to £{P — p, ±) also the third longitudinal polarization vector 



^(^-^'°)-i 



^" " .(l-.T)F+,-p.,- 



2(l-x)P- 
satisfying^ e(0) • e*(0) — —1, e(0) • e*(±) = 0, and [P — p) ■ e{Q) = 0, the corresponding additional LCWFs are 



(33) 



V'lo (^' Pt) = — ^7^ -TT H^ {Lz = 0) , (34) 

V2(l - x)Ma 

^U^,P.) - (- + ^^H;; + -^^) ^/. iL. . +1), (35) 

v2Aia 

i'+oi^^PT) = [V'-oC^^'Pt)]* (t)/x {L^ = -1), (36) 

V'lol^^^Pr) = --0|o(2;,Pr) 0/2: (Lz = 0). (37) 

From the above combinations we deduce, for example, that the proton with positive helicity +i can be in a state 
with probability density proportional to |'0-Pj where the quark has opposite helicity and L^ = +1 orbital angular 
momentum with respect to a scalar diquark. This configuration is relativistically enhanced with respect to |'0+ P with 
Lz = 0, where proton and quark helicities are aligned; thus, it suggests a possible explanation of the proton "spin 
puzzle" in terms of the relativistic aspects of the motion of quarks inside hadrons |5l| . 

For the purpose of this work, it is also important to note that a nonvanishing relative orbital angular momentum 
between the quark and the diquark implies that the partons do not necessarily occupy the lowest-energy available 
orbital (with quantum numbers J^ — ^ and L^ — 0). Hence, in this version of the spectator diquark model 
the nucleon wave function does not show a SU(4)=SU(2)(X)SU(2) spin-isospin symmetry, contrary to what is usually 
assumed [42]. 

Finally, we mention that the completeness relation for the last choice of the polarization sum in Eq. pH)) should be 
written 

Y, e*^{P~p,\a)e^{P~p,Xa)'e*^{P~p,t)e^{P~p,t)^-g>^'^, (38) 

where the unphysical time-like polarization state e^{P — p,t) = {P ~ p)^^ /Ma appears. The associated LCWFs read 
,4-, X p'i + xM^ — mM (1 — x)'^ ,, /^ „x /„„x 

i^U^,PT)^^ — °^ , ,;^ '-Hx (i. = 0), (39) 

V2(l + X)Ma 

^pZt{x,Pr) = y=— '-(p/x (L, =+1), (40) 

v2Ma 

i^:^t{x,pT) - [■tljlt{x,pr)]* (p/x {L, = -1), (41) 

ipZt{x,pT)^-'4}Xtix^PT)4>/x (iz=0). (42) 



Note that e(0) is not parallel to (P — p) because it describes longitudinal polarization states in the light-cone. 



C. T-even functions 

The simplest example of T-even parton density is the unpolarizcd quark distribution fi(x,pT), defined as 
/i(i,Pr) = J Tr [(»(i,p„ S) + <|i(i,p,, -5)) 7+] + h.c. 



h.c. 



By inserting in Ai ^°-' of Eq. ([9|) the rules (fTTjl for the nucleon-quark-diquark vertex, the dipolar form factor of Eq. (J13p , 
and the first choice in Eq. (jlOp for the sum of the polarization states of the diquark (transverse polarizations only), 
we get 

Ma), ^^ 9I [pI (1 + ^') + {m + xMf (1 - xf] (1 ~ x) 

h i-'P-)-(2,)3 2[pl + Ll{Kl)Y ■ ^^'> 

The same result can be recovered through the alternative definition 

Aiv— i Ag — ± 



Ajv = ± A„ = ± A„=± 



and replacing the results for the LCWFs using Eqs. ((20|) and (|32p for the scalar and vector diquark, respectively. 

If we use, instead, the second option of Eq. P^ for the sum over polarizations of the vector diquark (transverse 
and longitudinal polarizations), we obtain 

/f('^'(x,p,) + ^(|V^+oP + lV^+oP)- (48) 

The complete expression is given in Eq. (|A22ll and corresponds to Eq. (10) of Ref. [36] with Rg = 0. 

Finally, the results with the last choice of Eq. ([TU]) (transverse, longitudinal, and time-like polarizations) can be 
written as 

ff^\x,p.) + Yi(l^|oP + l^+oP) - I^(lV^ltP + \^^-t\') ■ (49) 

Note that the contribution of the diquark time-like polarization states enters with an overall negative sign. The 
complete expression is given in Eq. (jA26p and corresponds to Eq. (8) of Ref. [4q . 

Turning back to our preferred choice, i.e. the first option of Eq. (jlOp (light-cone transverse polarizations only), we 
now compute all other T-even, leading- twist TMDs. Their definition in terms of traces of the quark-quark correlator 
can be derived from, e.g., Eqs. (3.19) and ff. in Ref. [15]. To write them in terms of LCWFs, we need to introduce 
the polarization state in a generic direction St = {cos (jjs, sin (jjs) in the transverse plane 

U{P, T) = -^ (C/(P, +) + e*^-^C/(F, -)) , (50) 

C/(F, i) = -^ (C/(P, +) + e''(^''+-)c/(P, -)) . (51) 

For (pg — 0, 7r/2, we recover the (positive) polarizations along the x and y axis, respectively [53j. For the quark, we 
will use similar decompositions and use the notation Sqt and (j)s , i.e., 

uip, T) = -^ {u{p, +) + e~''^'^u{p, -)) , (52) 

u{p, i) = -^ (^(P' +) + e-*(^-''+")^(p, -)) . (53) 
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With this conventions and keeping in mind that Ax is absent for the scalar diquark and Xx — ± for the vector 
diquark, we can write the TMDs in the following way 



9lL( 



(-'P-)-T7^E(|V4aJ^-|V' 



-Ax I 



Pt ■ "Jt 

M 

W 



Xx 
Ax 



T |2 L/,T |2 

+Axl 



\^P[ 



+ |2 _ L/,+ |2 



\^ 



iA^ 



St ■ SqT hiT{x,pT) + -^-jj-^ "" j,j '^^ h^j,{x,pT) = 7— J Yl \\^ 



M M 



Ax 



T |2_L/,T |2 



IV; 



iA^ 



The above results automatically fulfill positivity bounds Q. 
The explicit expressions are 



9lL (x^Pt) 






9(s) 
9lT 



t/lT 






/.^ 



q(a) 



'it 






'it 



.g2 [(m + xM)2 - pI] (1 - a;)^ 



x,pI) 



(27r)3 2[p2+L2(A2)]4 

52 [p2 (1 + ^2) _ (^ ^ ^^,j)2 (1 „ ^)2] (1 „ ^) 



(27r)3 



2[p2+L2(A2)]4 



a;,Pr) 



52 M {m + xM) (1 - a;)3 

(2^ [p2+L2(A2)]4 



x,p^) 



a;,p^) 



a;,Pr) 



2 , _ 5^ xM (m + xM) (1 - x)2 



(27r)3 [p2+i2(A2)]. 



gl M (m + xM) (1 - a;)^ 



(2^)3 [pl + Ll{KW 
gl M (to + xM) (1 - a;)^ 



(27r)3 [pl + Llikl)]^ 



x,pI) = 



x,pI) 



gl [pI + {m + xMf] (1 - xf 



(27r) 



9l 



2[pi + mM)Y 

p^.x(l - x) 



(2^)3 [p2+i2(A2)]4 



a;,p^) 
x,pI) = 



gl AP (1 ~ xf 



(2^)3 [p2+i2(A2)]4' 



From the last two formulae we deduce also the expressions for the transversity distribution: 



hf\x,pl) = hlP{x,pl) + ^ ht^'^'\x,pl) 



gl (to + xMf (1 - a;)3 



(27r)3 2[p2+L2(A2)]^ 



hf^\x,pl) 



gl pI x{1 - x) 



(27r)3 [p2+L2(A2)]. 



(54) 
(55) 
(56) 
(57) 



(58) 
(59) 

(60) 
(61) 

(62) 
(63) 

(64) 
(65) 

(66) 
(67) 

(68) 
(69) 



Note that the functions gix and hj^j^ arise from the interference of LCWFs with \Lz\ = 1 and L^ = 0. The function hj^j, 
requires the interference of two LCWFs that differ by two units of L^ ■ This condition is necessary but not sufficient 
to have hj^rp ^ 0. In fact, the vector diquark spectator gives hj^j. = even if LCWFs with L^ = ±1 are present. 



Some interesting relations can be evinced from the above expressions. For example, the transversity with scalar 
diquark saturates the SofFer bound, while for axial-vector diquarks the relation is more involved: 

hf\x,pl) = I {ff'Hx,pl)+gf\x,pl)) , (70) 

hf^\x,pl)^-^l {ff''\x,pl)+gf''\x,pl)) . (71) 

When restricting to the results with scalar diquark, the giT distribution is connected to two other partners by the 
relations 

gt^\x,pl) = -ht^^'\x,pl) , g'Pi^^Pl) -^^hf\x,pl) , (72) 

171 + xM 

while for axial- vector diquarks we have 

gf^\x,pl)=xh^^^^\x,pl). (73) 

This relation is however different when considering spectator diquarks with more degrees of freedom (see App. VK\ . 
Our results seem to indicate that no general relation exists between giT and ^j*^^, contrary to what is proposed in 
Rcf. 4^. The reason is connected to the difference between LCWFs with Lz = \ and Lz — —1, as in Eqs. p4|) and 
([25|l . We also observe that in the vector-diquark case giL — hi and h^j. are not simply related through the relation 
suggested in Ref. [SJ] ■ We are led to conclude that such a relation is not general. 
The Pr-integrated results are 

Ms) ( .^L [nm + xMr + LliAl)]{l~xf 

^1 ^^"(27r)2 24L6(A2) ^'^> 

Ma),, gl [2 (m + xMf (1 - xf + (1 + x") Ll{t.l)\ (1 - x) 

^' ^^"(27r)2 24Lg(A2) ^'^^ 

„.Wm~ 5' {2{jn + xMf-Ll{t.l)\{\-xf 

,(a). .__ al [2 (m + xMf (1 - xf - (1 + x^) £2(Ag)] (1 - x) 
'' ^^" (2^)^ 24L6(A2) ^''^ 

hi(s)i._ gl {m + xMfjl-xf 
' ^^"(27r)2 12L6(A2) (78) 

' ^^" (2^)^ 12L^.(A^) • ^''^ 

D. T-odd functions 

The two leading-twist T-odd structures are the Sivers and Boer-Mulders distributions. They are defined as 

'JP^ /ii(x,p2) = _ 1 Tr [mx,p^, S) - <i>{x,p^,-S)) 7+] + h.c. , (80) 

^^ hi{x,pl) = - Tr [mx,p^, S) + ^x,pr, -S)) za'+ 75] + h.c. . (81) 

At tree level, these expressions vanish because there is no residual interaction between the active quark and the spec- 
tators; equivalently, there is no interference between two competing channels producing the complex amplitude whose 
imaginary part gives the T-odd contribution. We can generate such structures by considering the interference between 
the tree-level scattering amplitude and the single-gluon-exchange scattering amplitude in eikonal approximation, as 
shown in Fig. [3] (the Hermitean conjugate partner must also be considered). This corresponds just to the leading- twist 
one-gluon-exchange approximation of the gauge link operator of Eq. ^ [5^ . 
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p- p 



FIG. 3: Interference between the one-gluon exchange diagram in eikonal approximation and the tree level diagram in the 
spectator model. The Hermitean conjugate diagram is not shown. 



For the moment, we use Abelian gluons. The QCD color structure will be recovered at the end. The Feynman rules 
to be used for the eikonal vertex and propagator are [ll, [Sa] 



P 



i-l) 



-1+ 



(82) 



where Cc is the color charge of the quark and the sign of ie for the eikonal line corresponds to the gauge link of SIDIS. 
In cut diagrams one must take the complex conjugate of these expressions for vertices and propagators on the right 
of the final-state cut. 

The explicit form of the contribution <I>(^^ to the correlation function corresponding to Fig. [3] is 



$«(x,p,,5) 



1 



1 fir.io) 



(27r)3 2{l-x)P+ 
where T(a;,Py) is defined in Eq. ([8]) and 



(^M^°\S) M^^\S) + M^^\S) M^°\S) 



P^=t{x,p } 



(83) 



M'^^\S) 



dH ie,V,pn'_{jl>- l + m)ysU{P,S) 



(27r)4 {Di + ie) {D2 - ie) {D3 + ie) {D4 + ie) 
dH ie,e%[P - p, Xa) r^>^(j^ - / + m) d^,[p -l-P)y^ U{P, S) 



f27r) 



(Di + ie) {D2 - ie) [Ds + ie) (L»4 + ie) 



scalar diquark, 



axial- vector diquark, 



where for convenience we have introduced the notation 

Di=P - 



9' 



D2 = l+, 

D4^iP-p + l)^- 



(84) 



(85) 



M 



X 



In order to explicitly calculate M^^\ we need to model the gluon vertex with the scalar (Fs) and axial vector (Fa) 
diquark in Fig. [3l 



TP = ieci2P-2p + l)P 

-ie, [{2P ~2p + l)p /" -{P-p+{l + Ka)l)^9; -{P-P- kJ)" 9; 



a p 



(86) 



where e, is the diquark color charge, which is the same for scalar and vector ones and identical to that of the quark; Ka 
is the diquark anomalous chromomagnetic moment. The structure of the vector diquark-gluon vertex resembles the 
one for the coupling between the photon and a spin-1 particle (see, e.g., Ref. [531); for k^ = 1 the standard point-like 
photon- IF coupling is recovered (see, e.g., Ref. |58jV 

The Sivers and Boer-Mulders functions can then be computed as 



C^T Pti^' 



Tj r±_ 



M 



flTix,PT)^- 



1 1 



1 



4 (27r)3 2(l-a;)P+ 



Tr 



'^■^''^^■hi{x,pl) = l ^ 



1 



M 



4 (27r)3 2(1 - x)P+ 



Tr 



(^M'^^'>iS)M^°\s) ~ M^^\~S)M^°\-S)') 7 
(a^(i) {S)M^"^ (S) + X^i) {-S)M^"^ (-5)) ^a'+ 75 



h.c. 



(87) 



h.c. 
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Again, results have been produced for the three different choices of both Eq. (|12p for the form factors at the 
nucleon-quark-diquark vertex, as well as of the axial-vector diquark propagator on each side of the diquark-gluon 
vertex in Fig. [3l Consistently with the case of T-even parton densities, here we show the results for the dipolar form 
factor of Eq. P^ and for the light-cone transverse polarizations of the vector diquark, i.e. the first choice in Eq. (|10p . 
the other combinations being listed in App. |B] Combining the rules (fTTj) with the ((86)) ones, we can rewrite Eq. ([87]) 
and (EHl) as 



P-Lg(s)^ 2^_ ffs 1 Mel {^-x? 

4 (27r)3 2{l-x)P+ [p2+L2(A2)]2 

1 Mel {l-xf 

4 (2^4(1-:e)P+ [pi + Ll{Kl)f 



ftT'''\^.pl) = -^ TTTZT^ .,: Ir.^ TZT^^WtT^ 2Im Jf (89) 

,i^''^^(x,p^) = ^^;^^7^^:i^f^^::^^^^^:^^^^2imJi\ (90) 



h'[''^^\x,pl)^ftT'^'\x,pl) (91) 

ht''^''\x,pl)^-lft/''\x,pl). (92) 

Note that for scalar diquarks the spectator model gives the same result for the Sivcrs and the Boer-Mulders functions, 
independent of the choice of the nucleon-quark-diquark form factor (see App. IBJ) . 

In Eqs. (|89p and ff., the expressions Ji contain the integral over the loop momentum, the denominators -Di,2.3,4 
defined in Eq. ([55)l . and the evaluation of the trace of the projected amplitude. For instance (see App. IB ip 



j« ^ / J-L 9s{{p If) ^ / + _ +\ / + _ + It^\ ^ ^ 

1 J (27r)4 {Di + ie) {D2 - is) {Da + ie){Di + ie) \ ^ ' J\ ^ pI ^ 

To calculate its imaginary part, it is sufficient to make the replacements 

^ 27riS{D2), ^ V ^ -27riS{Di) (94) 



D2 — ie D4 + ie 

which corresponds to applying the Cutkosky rules }59i] . cutting the diquark propagator {D^) and the eikonalized quark 
one {02). We then get 

^^"^^^ ^Ii^^ '^%dT ^ (^+ + 2(1 -.)P+) [l+M-P+im + xM]^-^^ {2m)S{D2){-2m)5{D,) ^^^^ 

= -4P+ (m + xM) (1 - x) gsli . 

The calculation of Xi depends on the form factor used. Their calculation can be found in App.lCl For the case of the 
dipolar form factor we obtain 

Ap+f , nr\(A \ Tdip P+ {in + xM) (1 - xf 

^^P^{m + xM)il^x)9.I,^ = g.-^^^j^^j^^^^^^. (96) 

If the T-odd structures were deduced from the Drell-Yan amplitude, the Ai'-^^ of Eq. ([M)l would involve a (Z+ -|- ie) 
propagator, leading to the opposite sign in the cutting rule for D2. In the spectator model, this is the origin of the 
predicted sign change for f^ and h^ when extracting them in Drell-Yan spin asymmetries rather than in SIDIS 
ones [16i|. Analogously to Eq. ([55)) . we obtain 

„ „4- / , ^N ^dir, 2P~^ x(l — x) (m + xM) 

2l^Jt^-,P-xim^xM)g.lt-^g. ,Ll{[inpi^LliKl)l ' ^''^ 

By inserting these results in the model expressions of Eqs. (|55|) to ([5^ . we come to the final form of the Sivers and 
Boer-Mulders functions with scalar and axial vector diquarks: 
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JlT K-^^Pt) — A /r,-\d r2/A2\r„2 , r2^A2M3 ^^^> 



Pi9(.)^ 2^_ ffs ^eg (l-x)3(m + xM) 

4 (2^)4L2(A2)[p2+L2(A2)p 

4 (2^)4L2(A2)[p2+^2(A2)]a 



/.^^^^^(a:,p2,) = AV^^^(:z:,p2,) (100) 

/i^«('^)(2;,p2j=-iAV''^)(:r,p2). (iQi) 



X 



To connect the "Abelian" version of the gluon interaction to the QCD color interaction we shall apply the replace- 
ment Tj] 

el -> A-KCpas. (102) 

The Sivers and Boer-Mulders functions obtained in our model behave as 1/p^ at high p^, similarly to the /i in 
Eq. (U?]). As observed also in Ref. [60], this leads to a breaking of the positivity bounds for sufficiently high 
values of p^. This problem is due to the fact that the T-odd functions have been calculated at order a^, while the 
T-even functions at order a'g. At high p^., QCD radiative corrections generate a 1/p^ tail for /i and a 1/p^ tail for 
fvr UM ■ Our model is supposed to be valid for p^ ~ AP and for reasonable choices of the parameters no problems 
with positivity occurr in this region. 

Often the following transverse-momentum moments of the Sivers and Boer-Mulders functions are used: 

ftT^'\^)^ldp^^fM^,pl) 
ftr^'^'\x)=Jdp,\^fM^,pl). (103) 

In our model, they turn out to be 



^1^ ^'~ 32 (2^)3 M [L2(A2)]2 ^'"^^ 

f^g(a)(i)^^N_ gg Sc x{m + xM){l~xf 



;^^9(.)(1)(^)^J^-L^<7(«)(1)(^) (106) 

;,^«(a)(l)(^) ^_l^±^9(a)(l)(^^^ (107) 



f±q{s){i/2)i ^ gl el {m + xM){\-xf 

256 (2^)2 [L2(A2)]V2 






,i^q{a){i/2),._ 9l el x{m + xM){l-xf 

^1^ ^ ^ " 256 (27r)2 [i2(A2)]5/2 ' ^'^^) 

^^.(^)(l/2)(^)^^X^.(.)(l/2)(^) (110) 

;,^.(a)(l/2)(^) ^_i^±^,(a)(l/2)(^)^ (HI) 

E. T-odd functions: overlap representation 

As already mentioned above, T-odd leading-twist parton distributions arise from the interference of two channels 
leading to the same final state; for the case considered here (and depicted in Fig. [31), the two channels are given by 
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the tree-level and the single-gluon-exchange scattering amplitudes, respectively. In Ref. [6l|, it was suggested that 
T-odd parton densities can also be represented by overlaps of LCWFs, as for their T-even partners, provided that a 
suitable operator is included to describe the final-state interactions (FSI) produced by the gluon rescattering. So far, 
this representation was fully developed in a spectator model only for the Sivers function with scalar diquarks [2J, |23| . 
Here, we generalize it to the case of axial-vector diquarks, as well as to the Boer-Mulders function. In this way, all 
leading-twist (T-even and T-odd) parton densities can be given by overlaps of LCWFs consistently within the model, 
contrary to the statement of Ref. [23| . 

Following Ref. [26], for a nucleon transverse polarization state described by Eqs. ([50]) and (|?T|) along a generic 
direction St — {cos 4>s, simps), and for an analogous quark state described by Eqs. ([5^ and ([55]) along SqT — 
(cos c, sin 05^), we can rewrite the Sivers (|80|) and Boer-Mulders (|8T|) functions according to the Trento Conven- 
tions [621 (keeping in mind that Xx is absent for the scalar diquark and Xx — ± for the vector diquark) as 



2 {St XPt)-P 

M 



Iit{x,pI) 



{SqT -x-Pt)- P 
M 



dpT 



G{x,Pt,Pt) 



Yl [^IrAx(-^'P^)V^I,Ax(^'PT) - V'i*Ax(2^'P^)^i,Ax(2^'PT) 



Aq.A^ 



hi{x,pl) = 



dpT 



G{x,Pt,Pt) 



\ E K"x (^' Pt) ^Ux (^' p't) - V-xaV (^' P-) ^'ix. (^' Pt)^ 



Ajv, A> 



h.c. 



(112) 



(113) 



The above equations should be considered as assumptions, since it is not known a priori if the FSI operator 
G{x,Pt,p't) can be isolated and is the same for all functions and all types of diquarks. In our model, it turns out to 
be actually the same in all cases. In order to determine it, we must insert here above the expressions for the LCWFs of 
Sec. Ill B[ and compare the results with the ones from Eqs. ((89)) to (jM]) . after replacing ImJf , IraJi with Eqs. (|95ll97p . 
respectively, while keeping the definition of Xj^'^ (see App.[C|). For the scalar diquark case, for example, we get 



J IT 



{^^Pl) = {h 



M (1 - xf (m -I- xM) 



[Pi 



rn^'s)? 



dpT 



linG{x,pT,PT) {Pt - Pt) ■ Pt 

[p;,2+L2(A2)]2 p2 



(114) 



The above expression is identical to Eq. (|89p . after inserting Eq. ([M]) and the definition (jCSp of X^*^ (with the harmless 
substitution Z^ ^-> —I't), provided that 



lmG{x,pT,p'T) 



2(2^)2 {pt-p'tY 



CpOis 



27r {pt - p't 



(115) 



in agreement with the expression of Ref. _26J|. Following similar steps, we recover the same result (|115l) also for the 
Sivers function with axial- vector diquarks, and for the Boer-Mulders function as well. The FSI operator G{x,Pt,Pt) 
is indeed universal and describes a rescattering via one gluon-exchange, which corresponds to the expansion at first 
order of the gauge link operator of Eq. (O. Note that in other versions of the model (see App. |B]) the FSI cannot be 
as simple as Eq. (IllSp , since we observe also a dependence on the vector diquark anomalous chromomagnetic moment 
Kq, which is absent in the above equation. This does not imply that the FSI operator is not universal, but simply 
that it could have additional parts that are not interacting with scalar diquarks and transversely polarized vector 
diquarks. 

We close this section by observing that in our model we can generalize the relation between the first pT-moment of 
the Sivers function and the nucleon anomalous magnetic moment k, suggested in Ref. [26j in the simple scalar diquark 
picture. In fact, we define k in terms of the Dirac form factor using the overlap representation for the nucleon matrix 
element of the spin-fiip electromagnetic current operator 5l| : 



« = ^F.(0) 



-^ E -"/^*r(-,p;,A„)*,-(.,p.,A„ 



k,n.Xn 



(116) 



P'=PT 



where the sum runs upon the number of Fock states k, the number of constituents n in each state fc, and their 
helicities A„. Since in the diquark model of the nucleon initially at rest {Pt = 0) there is only one Fock state with 
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two constituents, and the kinematics of the diquark is constrained to the one of the valence quark, the wave functions 
^ reduce to the usual LCWF 26]. The momentum conservation for the struck quark reads p'^ = pj. + {1 — x) q^- 
Distinguishing between k^'^^ and k'^^"' for scalar and axial-vector diquarks, respectively, Eq. (|116|) becomes 

^qis) ^ i f dPrdx 






-2 ri 1 - x (1 - xf (m + xM) 



d:,i-±y±—!±^!lLZ^^ / dxK''^^\x 



^lia) 



27r)2 J, 12 [L2(A2)]3 



(2^)2 7o 12 [L2(A2)]3 

By comparison with the first pr-moment of fu? and fij? in Eqs. ([55)1 and ([M|) . respectively, 

ftT'^'\x)^Jdp^f^/'\x,pl) 

_ gl MCpas (1 - xf (m + xM) 
g2 MCpas x{l- x)"^ (m + xM) 



(27r)2 [2L2(A2)] 



3 



(118) 



we deduce the relation 






/iV(2^)=-:^^Ci=^«.T^, (119) 



valid for both types of diquarks, from which we have 

1 o 

dx (1 - x) f^j?{x) = -- MCpas n" (120) 

2 

that generalize the findings of Ref. [2J, [2^] ■ 

III. NUMERICAL RESULTS AND COMPARISON WITH VARIOUS PARAMETRIZATIONS 

In this section, after fixing the parameters of the model by fitting some known distribution functions, we show the 
numerical results of our model for a few selected TMDs. 

A. Choice of model parameters 

In order to fix the parameters of the model, we try to reproduce the parametrizations of parton distribution functions 
extracted from experimental data. When doing this, however, we have to face the problem of choosing a scale Q^ at 
which our model can be compared to the parametrization. In principle, this scale should be considered as a further 
parameter of the model. However, we checked that the lowest possible value of Q^ is always preferred by the fit. This 
is not surprising, since probably the model is applicable to a very low scale, beyond the limit of applicability of the 
perturbative QCD evolution equations. Therefore, we have decided to compare it to a parametrization at the lowest 
possible value of Q^. 

For the unpolarized distribution functions /" and ff, we have chosen the parametrization of the ZEUS collabora- 
tion [63] (ZEUS2002) at Qq — 0.3 GeV^. This set of PDFs gives also an estimate of the errors, which is important to 
perform a x^ fit. Other parametrizations either do not reach such low Q^ or provide no error estimate. 
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For the helicity distributions 5" and gf, we chose the leading-order (LO) version from Ref. [6j (GRSV2000) at 
Q^ = 0.26 GeV^. Since this parametrization comes with no error estimate, we assigned a fixed relative error of 10% 
and 25% to the up and down quark distributions, respectively, which is reasonably similar to the error estimates of 
other parametrizations at higher Q^ (see, e.g., Ref. [65]). 

Finally, in order to perform the fit we arbitrarily chose to select from each parametrization 25 equally spaced points 
in the range x = 0.1 to 0.7. 

The free parameters of the model include the quark mass to, the nucleon-quark-diquark coupling gx, the diquark 
mass Mx, and the cutoff Ax in the nucleon-quark-diquark form factor, for X = s,a scalar and axial- vector diquarks. 
It turns out that in order to achieve a good fit we need also to make a distinction between the two isospin states 
of the vector diquark. Hence, we will use ga Ma, and A^, for the coupling, mass and cutoff of the vector isoscalar 
diquark with 1^ = (corresponding to the ud system), and g'^, M^, and A^, for the normalization, mass and cutoff of 
the vector isovector diquark with I3 — 1 (corresponding to the uu system) . 

In order to reduce the number of free parameters, we decided to fix the value of the constituent quark mass to 
TO. — 0.3 GeV. We checked that the results are not very sensitive to the value of this parameter. 

To perform the fit, we need to discuss the relation between the functions /f , ff and ff , computed in the 
model, and the functions /" and ff of the global fits. For ease of interpretation, it is better to use normalized versions 

of the fi . Therefore, we write fi^oi-m = {^x/9x) /i where Nx are normalization constants determined by 
imposing 



- dx dplff^2^{x,pl) = 1. 

Jo JQ 



(121) 



Quite generally, the relation between quark flavors and diquark types can be written as 



fu _ 2 Ms) 2 Ma) (^0')\ 

Jl ^ '-s J 1 norm ' '-a J 1 norm \^^^) 

fi=<'ff::L. (123) 

We will refer to the coefficients cx as "couplings" , although they differ from the original couplings gx by the normal- 
ization constants Nx- They are free parameters of the model. 

In past versions of the spectator diquark model [42] , the quarks were assumed to occupy the lowest-energy available 
orbital with positive parity {J^ = ^ ); in this case, the proton wave function assumes an SU(4)=SU(2)(g)SU(2) spin- 
isospin symmetry, leading to probabilistic weights 3:1:2 among the scalar isoscalar (quark u with diquark s), vector 
isoscalar (quark u with diquark a), and vector isovector (quark d with diquark a') configurations. Moreover, the overall 
size of the couplings was adjusted to give a total number of three quarks. These choices led to the relations [42[ 

f" — _ /'"('') 4- _ f"('^) C1241 

J 1 c\ J 1 norm ' c\ J 1 norm V / 

fi = /itol • (125) 

There are two reasons to criticize this choice. First of all, in the present work the quark-diquark system can have a 
nonvanishing relative orbital angular momentum, as shown in the previous section. Thus, the proton wave function 
no longer displays an SU(4) symmetry. Secondly, strictly speaking the SU(4) decomposition gives coefficients that 
are three times smaller then the ones in the above relation. This is because the total number of quarks "seen" in the 
spectator model is only one, since the other two are always hidden inside the diquark. This is actually a fundamental 
limitation of the spectator model, it is independent of the SU(4) choice, and in our opinion it has not been sufficiently 
stressed in the literature. The only possible way out is to consider the diquark not as an elementary particle, but as 
formed by two quarks that can be also probed by the photon (see, e.g., Ref. |6a|). 

A different way to see this problem is by considering the (longitudinal) momentum sum rule. Since also the diquarks 
can carry momentum, they should be included in the corresponding sum rule.^ Using the handbag diagram of Fig. [2l 
we calculated the corresponding diquark distribution function /^ for the active diquark in the state X and the 
spectator quark with flavor q, again using the first choice in Eq. (fTUl) (independently of the choice of form factor). We 
found the remarkable property 

j'^'^\x)^!t''\\-x). (126) 



^ A similar approach has been used in Ref. [25l | to verify in the spectator model the validity of the so-called Burkardt sum rule 67], which 
is related to transverse-momentum conservation. 
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By splitting the total proton momentum sum rule into the contributions of quarks, Pg, and of diquarks, Px, using 
the symmetry property (jl26p we get 



Pa 



Px = dxx 
'o 

r-l 
dx X 
/o 

r-l 



2 ru{s) ,. 2 Ma) , . ,2 rdia') , -. 



2 Mu) ,. 2 ra{u) r. ,2 ra' (d) / x 

^s JlnormV-^/ ^ '-■a JlnormV-^^ ' ^a JlnormV-^/ 



(127) 



^ dx 
Jo 



c^ r^"'^ (x) + c^ /■"(") (x) + c'2 f^^"'^ fx^ 



= cl + cl + c'^ . 

It is therefore impossible in our spectator model to fulfill at the same time the momentum sum rule and the quark 
number sum rule. 

Although from the fundamental point of view it is more important to satisfy the momentum sum rule, from the 
phenomenological point of view it is impossible to reproduce the parametrizations in a satisfactory way. We decided 
therefore to avoid imposing the momentum sum rule and let the fit choose the values of the parameters cx- 



f;(x) 




gr(x) 





FIG. 4: The distribution functions fi{x) (above) and g\{x) (below) for the up quark (left panel) and the down quark (right 
panel). Data are a selection of 25 equidistant points in 0.1 < a; < 0.7 from the parametrizations of Ref. [63|] (ZEUS2002) 
and Ref. i64| (GRSV2000) at LO, respectively (we assigned a constant relative error of 10% to g" and 25% to gf based on 
comparisons with similar fits [651). The curves represent the best fit (x^/d.o.f. — 3.88) obtained with our spectator model. The 
statistical uncertainty bands correspond to Ax^ ~ 1- 



In summary, we have 9 free parameters for the model. We fixjthem by fitting at the same time /{', ff at Q^ = 0.3 
GeV^ from Ref. 63], and g^, gf at Q^ = o.26 GeV^ from Ref. Q at LO. The fit was performed using the MINUIT 
program. A x^ /d.o.i. = 3.88 was reached. The results are shown in Fig. 21 In spite of the very high x^, the agreement 
is acceptable, except perhaps for the down quark helicity distribution. The error band is deduced from the covariance 
matrix given by MINUIT and represents the standard 1-cr statistical uncertainty (Ax^ = 1). The corresponding 
values for the various model parameters are listed in Tab. ID 
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Diquark 


Mx (GeV) 


Ax (GeV) 


ex 


Scalar s (uu) 


0.822 ± 0.053 


0.609 ± 0.038 


0.847 ± 0.111 


Axial-vector a (ud) 


1.492 ± 0.173 


0.716 ± 0.074 


1.061 ± 0.085 


Axial-vector a' (uu) 


0.890 ± 0.008 


0.376 ± 0.005 


0.880 ± 0.008 



TABLE I: Results for the model parameters with dipolar nucleon-quark-diquark form factor and light-cone transverse polariza- 
tions of the vector diquark: the diquark masses Mx, the cutoffs A.x in the form factors, and the cx couplings for X = s,a, a' 
scalar isoscalar, vector isoscalar, and vector isovector diquarks. The fit was performed using the MINUIT program on the 
parametrization of fi{x) from Ref. ;6^ (ZEUS2002), and of gi{x) from Ref. H (GRSV2000) at LO, reaching a x^/d.o.f. = 
3.88. 



fi" (X, p?) 




0.1 0.2 0.3 0.4 0.5 0.6 



Pt (GeV^) 
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FIG. 5: The p% dependence of the unpolarized distribution fi{x,p%) for up (left panel) and down quark (right panel). Different 
lines correspond to different values of x. The downturn of the function /" at relatively small x is due to wavefunctions with 
nonzero orbital angular momentum. 



B. Unpolarized parton densities 



With the above model parameters, the proton momentum fraction Pq carried by valence quarks, is 



Pn 



dx X 



2 Ms) , , 2 ru(a) , . ,2 rd{a') , s 

^s JlnormV-^y T^ ^a J 1 norm V-^^ ' '^a JlnormV-^^ 

( dxx [f'^{x)+ ffix)] « 0.584 ±0.010, 
Jo 



(128) 



which is consistent with the ZEUS result of 0.55 [63J |. 

While for ff only the vector-isovector diquark plays a role, for /" it turns out that the contributions from the scalar 
and vector diquark have about the same size. The vector diquark is always dominant at high x. However, we know 
that the model is not reliable in the limit x -^ 1. In fact, the behavior at high x does not follow the predictions of 
Ref. [68| , since our model does not correctly take into account the dominant dynamics in that region. 

We consider now the p^ dependence of the unpolarized distribution function obtained in our model. In Fig. [5] we 
show the behavior of the up and down components as functions of p^ for a few values of the variable x. 

First of all, we observe that /" displays a nonmonotonic behavior at a; < 0.02. This is due to the contribution 
from LCWFs with nonzero orbital angular momentum. Although the details of where and how this feature occurs 
is model-dependent, it is generally true that the contribution of LCWFs with one unit of orbital angular momentum 
falls linearly with p^ for p^ ^ 0. This behavior is sharply different from the contribution of LCWFs with no orbital 
angular momentum. This simple example shows how the study of the p~ dependence of unpolarized TMDs can 
therefore already expose some effects due to orbital angular momentum. 

Finally, we observe that in our model the average quark transverse momentum decreases as x increases, and that 
down quarks on average carry less transverse momentum than up quarks. Although this is just a model result, a 
general message can be derived: the widely used assumption of a flavor-independent quark transverse momentum 
distribution is already falsified in a relatively simple model (see also Ref. |69|). 
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C. Longitudinally polarized parton densities 



The model parameters of Tab. |T] produce the axial charge 

gA= I dx [g'lix) - gfix)] = 0.966 ± 0.038 , 
Jo 



(129) 



in excellent agreement with the value 0.969 ± 0.096 deduced from the GRSV parametrization [6J|. 

It is, however, evident from Fig.[4]that our description of the down quark helicity distribution is in bad disagreement 
with the GRSV parametrization at large x. Nevertheless, we point out that there is a qualitative agreement with the 
parametrization of the so-called BBS model of Ref. [70] and the analogous parametrization of Ref. 68]. In particular, 
our model shows the same feature highlighted in this latter reference, namely that the contribution of the LCWFs 
with nonvanishing orbital angular momentum is dominant at high x. This is true in all distribution functions, but 
becomes particularly evident for the down helicity distribution, since the contribution from the LCWFs V'++ ^-nd 4'+- 
(carrying nonzero orbital angular momentum) are positive and make the distribution positive at x > 0.5. 



up, x=0.02 



down, x=0.02 
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FIG. 6: The p% dependence of the distributions fi{x,p%) — gi{x,p%) (solid line) and fi{x,p%) + g\{x,p%) (dashed line) for up 
(left panel) and down quark (right panel), at x = 0.02. The difference in their behavior is due to the different role played in 
the two combinations by wavefunctions with nonzero orbital angular momentum. 

The effect of orbital angular momentum becomes even more evident when considering the p^ behavior of the helicity 
distribution function. As an illustration, we show in Fig. [H] the behavior of the combinations fi{x,p^) — gilxjp^) 
and fiixjp^) + gi{x,p%). In the case of the scalar diquark, LCWFs with one unit of orbital angular momentum 
are filtered by the first combination. In the case of the vector diquark, the situation is opposite. The down quark 
distribution is entirely given by the vector diquark, therefore the /i(a;,Py) -\- gi{x,p'^) sum clearly turns down to zero 
for p^ -^ 0. For the up quark, the situation is less clear due to the simultaneous presence of scalar and vector diquark 
contributions. However, at a; = 0.02 the vector diquark is responsible for the nontrivial shape of fi(x,p%) + gi{x,p'^). 

It is interesting also to investigate the p^ behavior of g" alone, shown in Fig. [T] There is a dramatic change of 
behavior for different values of x, due to the difference between the scalar and vector diquark components of the 
function. If the spectator is a scalar diquark, for p^ = 0, where the LCWFs with orbital angular momentum vanish, 
the spin of the up quark has to be parallel to that of the proton, thus gi [x, 0) > 0. At high transverse momentum, 
where LCWFs with L^ = 1 dominate, the spin of the up quark has to be antiparallel to that of the proton, thus 
gi {x, oo) < 0. The situation is exactly reversed in the case of the vector diquark. As is already visible in Eqs. ([55)1 
and (j59p . at high transverse momentum the vector diquark always dominates and gives a positive result. At low 
transverse momentum, the relative size of the functions iJc {^\ ) i^^ the denominator determines which contribution is 
dominant. At higher x the scalar diquark dominates and gives a positive 5" (a;, 0), while at lower x the vector diquark 
dominates and gives a negative 5" (a;, 0). 

Once again, apart from the details specific to our model, these examples show that the exploration of the p^ 
dependence of the unpolarized and helicity distribution functions can expose very interesting features of the inner 
structure of the nucleon, related in particular to orbital angular momentum. 
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FIG. 7: The p% dependence of the heUcity distribution gi{x,p%). Different lines correspond to different values of x. 
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FIG. 8: The transversity distribution xhi{x) for up (left panel) and down quark (right panel). Dashed (solid) line for the model 
result before (after) the evolution at LO using the code of Ref. [7l| up to the scale of the parametrization from Ref. [Ta ]. whose 
uncertainty band due to errors in the fit parameter is represented by the shaded area. 

D. Transversity 

In Fig. [SI the predictions of the spectator diquark model for the transversity distribution are compared with the 
only available parametrization of Ref. [72]. In the left panel, x/i"(a;) is shown, whereas xhf{x) is shown in the right 
panel. All the model results at the assumed original scale Qq = 0.3 GeV^ are represented by the dashed line. The 
solid line indicates the result after applying the DGLAP evolution at LO up to the scale Q^ = 2.5 GeV^ using the code 
from Ref. [71|] . The latter scale pertains the parametrization of Ref. [7^ , whose errors in the fit parameters produce 
the uncertainty band represented by the shaded area. The model is in reasonable agreement with the parametrization, 
with the maxima in the correct position and a somewhat too small result for the up quark at small x. It should also 
be kept in mind that the present data reach at most x « 0.4 [73, [TJ] and, moreover, the ansatz of Ref. [74I does not 
allow for a sign change.^ 

Interestingly, for the up quark the model predicts a change of sign at a; ~ 0.5. To our knowledge, no other model of 
transversity displays this feature (see Ref. [76] and references therein; see also recent calculations in Refs. |66l. pTTI. pTa] ) . 



We point out that new fits of the transversity distribution functions have been presented at some conferences |7d| but not published yet. 
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The reason for this sign change is that the contribution of the vector diquark is negative, as evident from Eq. ((79|) . In 
our model, at moderate x the scalar diquark contribution is dominant, whereas at sufficiently high x the contribution 
of the vector diquark becomes in absolute size bigger, thus leading to the sign change. Other versions of the diquark 
model, even with vector diquarks, may not show this property. This is already evident from inspecting the results 
(listed in the appendices) for different choices of the diquark polarization sum. We don't think that our model 
calculation should be trusted more than others. Nevertheless, it might be interesting to contemplate the possibility 
of a sign change when choosing a form for the parametrization of the transversity function in "global fits." 
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FIG. 9; Same as in the previous figure, but for the Pt dependence of transversity at a:: = 0.1. 

In Fig. [9l the same comparison is performed as in the previous figure, but for the Pr dependence of the transversity 
at X = 0.1, as it is deduced from Eqs. (|68l69p . Again, there is a reasonable agreement between model predictions 
and parametrizations but for the trend of the result for the up quark at \pt\ > 0.3 GeV/c. However, we stress that 
the comparison may be affected by the different scale of the model results {Q^ — 0.3 GeV^) and the one at which 
the parametrizations are extracted {Q^ = 2.5 GeV^). The proper evolution of the TMDs has not been considered 
yet. It is interesting to point out that in our model hi{x,pl,) changes sign at pr ~ 0.5 GeV. This is due to the fact 
that the vector diquark contribution is always negative and dominant at high p^ . For the down quark, we note that 
hf{x,0) = 0, because the vector-diquark contribution to hi is entirely given by LCWFs with nonvanishing orbital 
angular momentum. 



XfiT (X) 



Xf,^^"^(X) 





FIG. 10: The first pT-moment xf-^^j, (x) of the Sivers function; left (right) panel for up (down) quark. Solid line for the 
results of the spectator diquark model. Darker shaded area for the uncertainty band due to the statistical error of the quark 
parametrizations from Ref. [Ty], lighter one from Ref. [80]. 
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E. Sivers function 



In Fig. [TOl the x/j^^ (x) moment of the Sivers function, predicted using Eqs. (|104[) and (|105[) . is given by the 
sohd hne and it is compared with two different parametrizations of the same observable. The darker shaded area 
represents the uncertainty due to the statistical errors in the parametrization of Ref. ^701 , while the lighter shaded 
area corresponds to the same for Ref. [SOI- Left panel refers to the up quark, right panel to down quark. First 
of all, we observe the agreement between the signs of the various flavor components, which also agree with the 
findings from calculations on the lattice |81l] . Also the maxima are reached at approximately the same x ~ 0.3 as the 
parametrizations. Instead, the "strength" of the asymmetry (related to the modulus of each moment) is too much 
weaker for the down quark, while it seems reasonable for the up one. Again, it must be stressed that no evolution 
was applied in the displayed model results. 
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FIG. 11: The model result for the spin density of unpolarized quarks in transversely polarized protons (see text for the precise 
definition) in Pt space at x = 0.1. Left panel for up quark, right panel for down quark. The circle with the arrow indicates 
the direction of the proton polarization. 



According to the Trento conventions |62l | , we define the spin density of unpolarized quarks with flavor q in trans- 
versely polarized protons as 



fq/pr{x,PT) = ft{x,pl) 



JIT \■^^ Pt) 



(Pxp^ 



M 



(130) 



In a SIDIS experiment, typically P is antialigned to the z axis that points in the direction of the momentum 
transfer q. Hence, if the proton polarization is chosen along the x axis, the spin density (jl30p shows an asymmetry in 
momentum space along the py direction, whose size is driven by the Sivers function. In Fig. 1111 we show fq/pr (0.1, Pr) 
for q = u (left panel) and q = d (right panel). Since the Sivers function for the up (down) quark is negative (positive), 
the density is deformed towards positive (negative) values ofpy. This feature is in agreement with the lattice results 
of Ref. [81| and with the signs of the anomalous magnetic moments k'' j61|] . 



Boer-Mulders function 



In Fig. [13 the xh^ ^ '{x) and xh^ ^ ' '{x) moments of the Boer-Mulders function, as deduced from Eqs. (|106ll07p 
and piOlllip . are displayed in the left and right panel, respectively. The solid lines correspond to the results for 
the up quark; dashed lines for the down quark. For the Boer-Mulders function, the only available parametrization 
appeared recently in [22[ , but the overall normalization depends on a parameter w that cannot be fixed with available 
experimental information. Our result agrees in sign and shape with that extraction. The absolute values of our 
functions correspond to w « 0.3. We remark that there is full agreement between the sign of the u and d components 
and the aforementioned lattice calculations [81| , as observed also in a different version of the spectator model [S^ and 
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FIG. 12: The xh^ (x) (left) and xh^ ' ' (x) (right) moments of the Boer-Mulders function. Solid and dashed lines for up 
and down quarks, respectively. 

in the bag model [43]. This agreement seems to be a general feature, as argued in Ref. [2^. ^ 

In Fig. 1131 we show, again at a: = 0.1, the spin density of transversely polarized quarks with flavor q in unpolarized 
protons, related to the Boer-Mulders effect by [62 | 



/gT/p(a:,PT) 



1 



fl{x,p^) - /ii ''(x,p^) — ^ 



(131) 



where now the quark polarization Sq points along x. Since the Boer-Mulders function is negative for both flavors 
(see Fi g. [T ^. the related spin density is always deformed towards positive Py, again in agreement with the lattice 
results fSll. 
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FIG. 13: The model result for the spin density of transversely polarized quarks in unpolarized protons (see text for the precise 
definition) in Pt space at a; = 0.1. Left panel for up quark, right panel for down quark. The arrow inside the circle indicates 
the direction of the quark polarization. 



^ A different result for the sign of the down quark Boer-Mulders function was obtained in Ref. [4a |. probably due to a mistake in that 
calculation (see App. [B4t . 
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IV. CONCLUSIONS 

We have presented a systematic calculation of all leading-twist parton distributions in the nucleon in a diquark 
spectator model. We have generated the relative phase necessary to produce T-odd structures by approximating the 
gauge link operator with a one gluon-exchange interaction. All results have been presented in analytic form and 
interpreted in terms of overlaps of light-cone wavefunctions. 

We tried to extend and improve the spectator model calculations presented in Refs. [3a, H, 13 by considering 
several choices of the axial- vector diquark polarization states and of the nucleon-quark-diquark form factor. We listed 
the analytic expressions for all possible choices in the appendices. We critically reconsidered some of the limits of the 
model and the choice of model parameters used in the past literature. In particular, we showed that the spectator 
diquark model is not able to reproduce both the quark number and momentum sum rule at the same time, because 
the diquark is considered as a charged parton, hence active in the sum rules. We argued that the proton wave function 
does not show the usual SU(4)=SU(2)®SU(2) symmetry [42|, since the quark-diquark system in its ground state can 
have a nonvanishing relative orbital angular momentum. 

For numerical studies, we chose the version of the model that in our opinion is more sensible and practical, i.e., the 
one where only light-cone transverse polarizations of the diquark are present and a dipolar form factor is used. We 
identified nine free parameters of the model and we fixed them by reproducing the phenomenological parametrization 
of unpolarized [63|] and longitudinally polarized [6J| parton distributions at the lowest available scale, i.e. Q^ = 0.3 
and 0.26 GeV^, respectively. 

Whenever possible, results have been compared with available parametrizations. For the chiral-odd transversity 
distribution, there is only one available from Ref. [73], which was deduced from SSA data at Q^ = 2.5 GeV^. The 
Pr-integrated model result, once evolved to this scale using the code from Ref. [7l(] at LO, displays a satisfactory 
overall agreement. The / , ^ (a;) moment of the chiral-even T-odd Sivers function /j^ was compared with the 
parametrizations of Refs. 79, 80]. There is agreement between the signs of the various flavor components and 
between the positions of the maxima in x, but the absolute value of the function is somewhat too small for the down 
quark. The comparison is affected by the difference of the scales, since evolution equations for the Sivers function 
have not been used. We also plotted the hi {x) and h-^ (x) moments of the chiral-odd T-odd Boer-Mulders 

function hj^ . We have also shown the quark spin densities defined in the Trento conventions [62l | , as produced in turn 
by the Sivers or the Boer-Mulders effects. For unpolarized quarks in transversely polarized protons, the spin density 
fq/pi is linked to f^, while for transversely polarized quarks in unpolarized protons the /„?/« is linked to hj^. For 
transverse polarizations along the x axis, the contour plot in the quark momentum space of such densities at x = 0.1 
displays a distortion in the Py direction, whose sign is consistent with the lattice findings for the corresponding spin 
densities in impact parameter space |81| . 

Using the model parton densities discussed above, various spin, beam, and azimuthal asymmetries in semi-inclusive 
hadronic reactions can be predicted, which are of interest for several experimental collaborations. Model calculations 
can be useful to interpret experimental measurements, helping us to explore spin-orbit parton correlations inside 
hadrons and shed light on the well-known puzzle of the proton spin sum rule. 
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APPENDIX A: T-EVEN FUNCTIONS IN DIFFERENT VARIATIONS OF THE MODEL 

In this Appendix we list the leading-twist T-even parton densities obtained in the context of our spectator diquark 
model, for all the choices of axial-vector diquark polarization sum and nucleon-quark-diquark vertex. To avoid 
overloading the notation, we will use the same ones for the parameters involved {gx, Mx, Ax). However, it must be 
kept in mind that the numerical value of these parameters can be different in the various versions of the model. 



24 
1. Scalar diquark 

The results for the scalar diquark are 

^1^ ^ '^"^" (27r)3 2[p2+L2(^2)]2 ' ^^^J 

.7^<^Vrn^l- -^^ M{m + xM){l-x) 

hf,^^'\x,pl) = -gf^\x,pl), (A4) 

/^?^^^(2:,p'J = /f^H:^,P^), (A5) 



(27r)3 [p2+i2(^2)]2 



"it (^;Pr) — ToZvT r„2 I r2('™2M2 ' (^o) 



where we recall that M is the nucleon mass and m is the mass of the parton. From the latter two densities, we 
construct the contribution of the scalar diquark to the transversity: 

hf\x,pl) = hf4\x,pl) + ^ ht^^^'\x,pl) 



(A7) 
__ g2 {m + xM)Hl-x) _ 1 / (,) , (,) , n 

" (27r)3 2[p2+L2(m2)]2 ^ 2 V' [^,Pt)+9i [x,p^) ^ 

The above results are valid for a point-like nucleon-quark-diquark coupling. For the other form factors it is sufficient 
to apply the replacements 

\p2 j^ ^2/^2\l2 

g^ -^ gK'^-xf I r2/A2M4 dipolar form factor, (A8) 

5,2->g^e-[p|+ix('"')]/[(i-=>:)Ax] exponential form factor. (A9) 

The integrated results are obviously different for the three form-factor choices. In all cases the transversity function 
saturates the Soffer bound, i.e., 

hf\x)^^[ff'\x)+gf\x)). (AlO) 

• Point-like coupling (to avoid divergences we assume that the p2 integration is extended up to a finite cutoff A2) 

,(s) _ 9l (1 - X) ("^ + ^Mf A2 - L2(m2) Ag + £2(^2) [A2 + L2(m2)] log (^ + l) 

•^1 ^^^ " (2^)2 4X2(^2) [A2+i,2(^2)] ' (^11) 

,(,) _ gl (1 - x) ("^ + ^^--^)' ^l + -^^("^') Ag - I^li^^) [A^ + ^2(^2)] log (^ + l) 

^1 ^^^ " '(2^)2 4X2(^2) [A2 + i2(m2)] ' ^^^^^ 

• Dipolar form factor [same as Eqs. ((74)) and ((76|) ] 

,,(s). . .9? [2(m + xM)2 + L2(A2)](l-x)3 

^1 ^"^^-(2^ 24L6(A2) ^^^^^ 

„.(^).,^ .9.' [2(m + xM)2-L2(A2)](l-.)3 

^1 ^"^^"^2^ 24Lf(A2) ^^^^^ 
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Exponential form factor 



ff'\-) 



gl 1 



l-x 



(27r)2 4 



2 Lf (™^)/[(l-x) A^] ^^ [(^ ^ ^^^)2 _ ^2(^2)] 



r 0, 



Ll[ni^) 
2X2(^2) \ 2 [(m + 2;M)2 - LUm')] - {\ - x) A^ 



(l-x)A2 
l~x 



A? 



5f ^W = (i)2 \ e--^(™^)/[(-^)A^l ^^ [(m + .xM)2 + L2K)] 



2L2(,„2) X 2[{m + xMf + Ll{m'^)]+Kl{l-x) 



where F is the incomplete F function, defined as 

/•OO 

r(a,z) = / f-^e'^dt 



A? 



(A15) 



(A16) 



(A17) 



2. Axial-vector diquark with light-cone transverse polarization only 

The unintegrated parton densities are 



gfT\x,pl) 



gl p^^ {I + x^) + [m + xMy {I - xY 



(27r)3 2[p2 +L2(^2)]2(l_2.) 

gl pI {I + x"^) - [m + xMf {I - xf 

(27r)3 2[p2 +L2(„j2)]2(l_2;) 

gl Mx{m + xM) 



(2^)3 [p2+L2(m2)]2 
h'^^'~^\x,pl)=gp{x,pl)lx, 



""IT V'^iPt/ 



,-L9(a) 



.9^ 



XPj 



(27r)3 [p2+i2(^2)]2(l_^) 



/i,^''^"^(:z:,p^) = 0, 



L9(a) 



<?(a)/ 



/.r(:z:,p^) = /ilV^(:.,p^). 



(A18) 



The above resuhs are vahd for a point-Uke nucleon-quark-diquark coupUng. For the other form factors it is sufficient 
to apply the replacements in Eqs. (|A8[) and (jA9[) . 
The integrated results are 

• Point-like coupling (to avoid divergences we assume that the p^. integration is extended up to a finite cutoff A2) 



ft''\x)-^ 



1 



(2^)2 4L2(^2)[A2+^2(^2)](l_^) 
+ L2(m2)[A2+L2(^2)j(^^^2)i^g 



a; A2 [(Af2 _ ^2) (1 - a;2) + 2mM (1 - xf - Ml (1 + x^)] 
A2 



Ll{m? 



gf^\,)=9^^ 



Li{m')[Ai + Liim')]il + x')log 



A2 



(2^)2 4L2(^2)[A2+i2(^2)](l_^) pv .La. ay uy ■ , ^\Ll{m^)' 

- A2 [(1 - x) mM [x{l - a;) {2M - to) + 2m] + x(l - x) M^ {x - 2x'^ - 1) + x{l + x^) MM^ 
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1 ^^>~ (2^)2 



A 



^[(1 - x) [xAP - m') - xMl] + Llim^) [A^ + LUm^)] log (^^ + 1 



Dipolar form factor [same as Eqs. (|75I77[ [75)1 ] 

gl [2 (m + :rM)2 (1 - x)^ + (1 + x^) ^^(a^)] (i _ ^) 



^1 ^"^^-(2^ 24L6(A2) ' 

,(a). gl [2{m + xM)^l-xf-{l + x')Ll{Al)]{l-x) 

^' ^ ' (27r)2 2ALUAI) 



i <2(a) 



/i?^"'(a;) 



.9a 2;(1 - x) 



(27r)2 12L4(A2) 



• Exponential form factor 



ff''\x) 



A^L^K)(l-x)(l + x^)r 0,2 



(27r)2 4A2L2(m2)(l-x)2 ^ 

+ [(m + xMf (1 - .t)2 - Llim'^) (1 + a;^)] 



2^^2^ 



Li{m 



{l-x)Al 



9t''\-) = 7^. 



(1 -x)Ale -2^^(™^)/[(i-) A^l - 2L^K) r (o, 2 (^^^^) 



AiLi{m')il-x)il + x')T[0,2 



il-x)Al 



(27r)2 4A2L2(m2)(l-.x)2 

+ [(to + xMf (1 - xf + Ll{m^) (1 + a;2)] 



/if"^(a:) 



9l 



(27r)2 2A2 (1 - xy 



A^ (1 - x) 



-2Ll{7n^)/Hl-x)Al] 



r 0,2 



(l-a;)A2 



2L^K)r 0,2 



{l~x)Al 



(A19) 



(A20) 



(A21) 



3. Axial-vector diquark including also longitudinal polarization 

The unintegrated parton densities are 



ff'\x,pl) 



gl 



1 



pi, + xM^ {2pl + xM^) 



{27rr A[pl + Ll{m^)]^MUl-x) 
+ (1 - xf [pI [A'f + rr? + 2M^) + 2w?MI + GxmMM^ + 2x^M^M^ + m^M^ (1 - xf] 



2 

Ma)( 2\ _ 9a 



9iri^,pi) = 



pI + xMI{2pI + xMI) 



?(") t 2 \ 



{2^Y^pl + Ll{ra^)YMl{l-x) 
+ (1 - xf [pI {2MI - TV? - M^ - AmM) - 2m^Ml - 2xmMMl - 2x'^M^M^ + m^M^ (1 - xf] 

gl M (to + M) pI - mM (m + M) (1 - xf + xM^ [M{2x - 1) + to] 



(2^)3 



2[p2+L2(TO2)]2M2 
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,-L9(a) 



{x,pI) 



glM (m + M) pi - mM (m + Af ) (1 - x)'^ + xMM^ - mMl{x - 2) 



"it I'^iPrJ 



(2^)3 2[pl + Ll[m^)fMl 


' 


gl pi + [2xMl + (m2 + A.p) [l - xf]pl + [xM^ + mM (1 - 


-x)r 



,-L9(a) 



{x,pI) 



gl M^ (to + Af )2 (1 - x) 



hf''\x,pl) = 



(2^)3 2[p2+L2(m2)]2M2 ' 

52 p4 ^ [2xA//2 - 2TOA-f (1 - a;)2] pi + [xM^ + mM (1 - a;)2]2 
~ (27r)3 4M2 (1 - x) [p2 + a;M2 + (1 - a;) (to2 - a;Af2)]2 



(A22) 



The above results are valid for a point-like nucleon-quark-diquark coupling. For the other form factors it is sufficient 
to apply the replacements in Eqs. (jASP and (jAOp . 
The integrated results are 

• Point-like coupling (to avoid divergences we assume that the pi integration is extended up to a finite cutoff' A2) 



f,^''\x) 



9l 



Kl[Ll{m^) + KlLl{m^)] 



(2^)2 %Ll{m^)[Kl + Ll{m^)]Ml{l-x) 
+ {l- xf [2Ml [m^ - Ll{m^)] + x [{m^ - M^f - M^ (to^ - 6mM + M^)] + x^ 2M^M^] 

+ Llim') [Al + Llim')] (1 - x) log ( -^ + l) [(1 + x) (m^ - M^) - 2Mf (1 - x)] 



L2(to2) 



5f-)(,)^ 5a 



1 



xMLiim') At [Ai + L^ („,/)] + A^ (1 - xY 



(2^)2 8i2(^2)[A2+i2(^2)]MAf2a;(l-x) ^ 

Ll{m^) [mM^ {2x - 1) + m^M {x - 2) - m^ + xM {M^ - 2M^)] 
-H TO^ (1 - a;) - TO-^Af (a; - 2) + m^M^x + m^M"^ (1 - x'^) + m^M^x {x^ - 2a; - 1) 
+ mM^M^x + mM'^x {2x'^ - a; - 1) + x^M^ (Af ^ - 2Af^) 

A2 



+ il-x)Li{m')[Ai + Li{m')]lo. 



4-1 



L2(to2) 

Ll{m^) (m + A/) - m^ (1 - x) + m^M {x^ -2x-\)~xm (SAf ^ (1 - a;) + Ml) 
-H xM {Ml(\ - 2x) + 2xM'^) 



ht\x)^--^^ 



1 



L2(to2)A^(L^(to2)+A1) 



(27r)2 8X2(^2) [A2 + Ll(m^)\ Ml (1 - x) 

-F (1 - xf A\ [(m + Mf Ll{m^) +x{m- Mf ((to + Mf - Ml)] 

/ A2 
- 2Ll{m^) {Llim^) -f A^) (to -|- M) (m - a;Af) (1 - x) log 



Llim^) 



(A23) 



• Dipolar form factor 



ff''\x) 



al 



\-x 



2Ll{Ai) + 2xMl{Li{Ai)+xMl) 



(27r)2 48L6(A2)Af2 
+ (1 - a;)2 [(2Af 2 + m^ + Af 2) i2(A2) + 4m2^^2 ^ ^ UmMMl -f a;^ 4A//2A//2 + 2m^M^ (1 - a;)^] 
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5f)(x) = 



2Ll{Ki) + 2xMt{Li{K)+^Mi) 



(27r)2 48L6(A2)M2 
+ (1 - xf [{2Ml - m^ - M^ - AmM) LI{AI) - Arr?Ml ~ x AmMM^ - x^ AM^M^ 



+2Tr?M^{l -xf 



hf''\x) 



9l 



l-x 



{2i.f2ALl{Kl)Ml 



[Ll{kl) + [xMl-mM{l-xf]Ll{kl) + [xMl+mM{l-xfY] . (A24) 



Exponential form factor 



lt''\^) 



9l 



1 



(27r)2 lQLl{m?)KlMl{l-xf 
X I 2 



V (1 - x) K 

(1 - xf \rr? {2Ml + M^ (1 - xf) - Ll{rr?) (2M^ + m^ + M^) + 2xMMl (xM + 3m) 



+ (1 - x) xM2 (xM2 - m') + Ll{m^) {Ll{m') - xM^) 
+ AlLlim')il-x) 



(1 - xf 2r ( 0, ,^^ 7a2 ) ("^^ + ^'^^ + 2^a) 



9l^''\-) 



9l 



{l-x)Al 
+ {l-x) Ale -2 [p'T+Llim-)]/[(i--) All + 4r f 0, ^^'("^'M [^Mf - Llim^)] 

1 



(27r)2 16L2(TO2)A2M2(l-a;)2 

A2 (1 - x) e -2 [Pl+il(™^)]/[(i-) All _ 2 l2(to2) r (o, -^^^^ 



(1 - x)"^ 2m^M^ - (1 - x)2 [2toM[xM^ - i^Cm^)] - Ll{m^) (m + M)^] 
+ (1 - x) [Al Llim^) + 2Ml {vn? (x - 2) + (2x - 1) xH'I^ 
+ 2Ll{m^) [Llim') - xMI - 2MI (1 - x)^] 

2Ll{m^) Al {l-x)T (^0, ^^^) [i^(™') - 2xMf + (1 - x)^ [{m + Mf + 2mM - 2M!]] \ , 



hl^''\x) 



9l 



1 



(27r)2 16L2(m2)A2M2(l-x)2 
2L2(m2) 



sinh(2^- + "-(-:)Vcoshr2^- + ""(-^) 



2LfK)r 0, 



"^^ U V (i-^)A2 ; V (i-2;)A2 

(1 - x)"* 2to2m2 + (1 - •^)^ 4toM [xM^ + ^^(m^ 



(l-x)A2 

+ (1 - x) [Al Ll{m') + 2xMl {xM^ - w?)] + 2Ll{m^) (Ll{m^) - xMf) 
2i^K)A rr2, 2 



+ 2Lf(m2)A2(l-x)r 0, 



(l-x)A2 



£2(^2) _ 2xM2 + (1 - x)2 2mM] 



(A25) 
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4. Axial-vector diquark including also time-like polarization 



The unintegrated parton densities are 



Ma), 2^_ gl [pI + {m + xM)^ + 2mMx] {1 ~ x) 
h [^.Pt) (2^)3 2[p2+L2(m2)]2 

l{a), 2^ ffg hpl + m^ + xHP] jl - x) 

9lL [^^Pt) (2^-)3 2 [p2 + L2(m2)]2 

„1M( 2x 5L ^^'^(^-^) 

9lT y^^PT) (2^)3 [p2+^2(^2)]2 ' 



(2^)3 [p2+i2(^2)]2' 

The above results are vahd for a point-Hke nuclcon-quark-diquark coupHng. For the other form factors it is sufficient 
to apply the replacements in Eqs. (jASP and (JA9|. The result for /i with dipolar form factor corresponds to that 
obtained in Ref. [4a |. 

The integrated results are 

• Point- like coupling (to avoid divergences we assume that the p2 integration is extended up to a finite cutoff A2) 
gl (1 - x) [(^ + Mf + 2mM - Mf] x Kl + hUm') [Kl + L2 (m2)] log (^^ + l) 



ff''\x) 



gf^'i^)-^ 



(2^)2 ALl{m^)[Kl + Ll{m^)] 

g2 (^ _ ^) Kl jLlim^) + m2 + a:2Af2) - ^2(^2) [A2 + L2(m2)] log (^^ + l) 

^W 4X2(^2) [A2+i2(^2)] 



1 ^ ^ (2^)2 2L2(^2)[A2+^2(^2)] ■ ^ ) 



• Dipolar form factor 



Ma)( ^ ^ ^L [^a(Aa) + 2 [(m + xM)^ + 2xmM]] (1 - xf 

^' ^''' (2^)2 24L6(A2) 

,(a), , 3^ [Lg(Ag)-2(m2+x2Af2)](l-x)3 

yi I2;j 



(2^)2 24L^(A2) 

.g2 raM X (1 — : 
"(2^ 6i6(A2) 

Exponential form factor 



'^f^(-)--A^^^^^?fe#- (A28) 



.9(a), . _ gl 1 I -2L\m^)/[{i-x)A''] [(m + xM)2 + 2mxM - £g(m2)] (1 - x) 
^1 ^^^-(27r)2 4\' L2(^2) 

- r f 2 . ^'^""'^ ^ ^ t^"" + "'^^^^ + ^'"''^^ " ^'^""'^1 - (1 - ^) A^ 



(1-^)A2; A2 

g,(a)(^)_ gg \^rfo,2 ^'^"^'^ ^ 2[m2+x2M2 + Lg(m2)] + (l-a;)A2 



(27r)2 4 V ' (1--^)A2; A2 
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APPENDIX B: T-ODD FUNCTIONS IN DIFFERENT VARIATIONS OF THE MODEL 

As a continuation of App. \^ here we list the Sivers and Boer-Mulders functions, namely the leading-twist T-odd 
parton densities obtained in the context of our spectator diquark model, again for all the combinations of diquark 
propagators and nucleon-quark-diquark vertices. 

1. Scalar diquark 

For scalar diquarks, we have 

r±q{s)f 2^ 9siP^) 1 Mel 2ImJf 

■^1^ ^^'^^' 4 (27r)3 2(1 - x)F+ p2 _ m2 

ht^^^\x,pl)^f^^'^'\x,pl), (Bl) 

where the Jf' integral is defined as 

9s{ip-l?) 



f crl 



(27r)4 {Di + ie) {D2 - ie) {D3 + ie) [D^ + ie) 
Tr[f 



- / + m) ( ^ + M) 75 ^ ( :^ + m) (2P - 2p + l)p n^ 7 

J (27r)4 [Di + le) {D2 - ie) {D3 + te) (D4 + te) \ ^ ' 

(l+Mei^PT^Sr, - P+{m + xM) e'^ It^St,) , 



with Di, D2, D^, D4 defined in Eq. (j85p . The imaginary part of J( can be extracted by using the Cutkosky cutting 
rules on the loop diagram of Fig. [31 which in the present case amount to put on shell the eikonalized virtual quark 
propagator D2 and the virtual scalar diquark propagator 1)4. The resulting S functions (see below) reduce the integral 
in Eq. (|B2p to a bidimensional integral in (Pl^- In general, for a n-dimensional integral / d"l lpf{l,p) the term Ip can 
be replaced by the expression pp{l ■ p)/p'^ . For the present case n = 2 and with the identification Ip ~ l^t, Pp — Pt,, 
we finally can write 

2ImJf= f ^iA^-fl 4 {1+ + 2(1- x)P+) (l+M-P+{m + xM)^-^) {2m) 6(02) {~2m) S{D,) 
J (2tt)'^ D1D3 ^ ' \ Pt J 

= -4P+ (m + xM) (1 - x) g,Ii . 

(B2) 

The explicit expression of Ii clearly depends on the choice of the nucleon-quark-scalar diquark vertex form factor. 

The Boer-Mulders calculation gives exactly the same results as the Sivers one, in the scalar diquark framework, 
since the relevant trace over Dirac-Lorentz structures is now given by 



Tr 



{i»-]i + m){If + M){i, + m){2P-2p+l)pnP_ia'+^5 ^ -'ii{l+ + 2{l-x)P+^{l+ M e^p^,-P+ {m + xM) e^ It^ ■ 

(B3) 
• Point-like coupling 

2 Im Jf =gj ^ -1- 4 (;+ + 2(1 - x)P+) (l+M - P+ (m + xM) ^-Zl^A {2m) 6(02) {-2m) 5{D^) 

J (277)" D1D3 ' ' \ Pi J 

AP+( ^ nr^n ^ tp-'- P+ (m + xM) (1 - x) , f Ll{m^) + pl \ 
= -4P {m + xM){l-x)gslf ^-.9. ^^ ^"^ L^Jm^) ) ' 

(B4) 
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where Xf ' ' is calculated in App. [C] 
Using Eq. ([8]), the final result is then 

f^«(^V. r,2^ ^ _5l MjI (m + xM) (1 - x) ( Ll{m')+pl 

h^'^'\x,pl)^ftr'^'\x,pl). 

Dipolar form factor. The final results, already given in Eq. (|98|) and (jlOOp . are 

gl Mel {I - xf {m + xM) 



/it {XiPt) ^ f 



4 (2^)4L2(A2)[p2+i2(A2)]3 

h^'''''\x,pl) = ftr'^'\x,pl) 
• Exponential form factor 



2lxnJl^g, 



(2^)4 Di D3 

= -AP+ {m + xM) (1 - x) gs II""^ 



4 (Z+ + 2(1 - x)P+) ( ?+Af - P+(m + xM) ^^^ 



-5s 



P+ (m + xM) (1 - x) 



■npl 



r 0, 



(1-.t)A2 



-r 0, 



(1-.t)A2 



and I^'^^ is calculated in App. [C] 
The final results are then 



(B5) 



(B6) 



{2TTi) 5{D2) {-2m) 5{Di) 



/: 



'^i('^)(^ „2^_ 9I Mel {rn + xM){l~x) ^_^^2^^Ll{ra^)]/\(i~x)Kl] 

IT K'^tPt) — . /o„^4 ^2 rr2/'™2\ 1 ^21 



4 (27r)4 p2 [7^2(^2) +p2 



r 0, 



(l-x)A2 



-r 



Ll{m^)+pl 
(l-x)A2 



/l^^(^Hx,p^)^/^-L/«)(^,p2) 



(B7) 



We have 



2. Axial-vector diquark with light-cone transverse polarization only 



/-/")(.,p^) = ^ 



ht'^''\x,pl) = ^''^'''' 



Mel 



2 Im Jf 

(27r)3 4(1 - a;)P+ p2 _ „^2 



Afe? 



2Imj; 



/(a) 



4 (27r)3 4(1 - x)P+ p^ -TV? ' 



where now the Jf and J(° integrals are defined as 



(4v.^T,)jr - / 



ga{{p-iY 



dH 
(27r)4 (Di + ie) (Da - ie) [D^ + ie) (154 + le) 



Tr 



(B8) 



- / + m) 7'^ 75 ( ^ + M) 75 ^7" 75 (^ + m) rf^^p - ^ - P) rf,a(P - p) 



and 



^tPtjJ '-'1 



[(2P - 2p + Op g"" -{P-p+il + ^aWg; -{P-p- nJY g^] n^ 7^ 

J^Z ga{{p-l?) 

(27r)4 (Di + is) {D2 - ie) {D3 + ie) {D4 + is) 



(B9) 



Tr 



- / + m) 7^ 75 ( ^ + M) 7" 75 (|( + m) rf^,(p -l-P) d,a [P - p) 



[{2P -2p + l)p g"^- - (P - p + (1 + KaWg^; -{P-p- kJY <] n^. ia'+^5 



(BIO) 
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and the explicit expressions of the d^,j{p — I — P) and da-a{P — p) structures are those expressed in the first Hne in 

Eq. unD. 

• Point- hkc couphng (to avoid divergences we assume that the p^ integration is extended up to a finite cutoff A^) 



2 fm J^ - -8P+X m + xM) ga if' = -ga ^ log °;,/ ,, ^ 



where Tf' ' is the same integral as in Eq. (jB4p but with the substitution Ls{m'^) ^^ La{m?'). 
Using again Eq. ([8]), the final result is 

.^q(a), 2.9l Mel x{m + xM) f Ll{m^) + pI 



4 {2nYpl[Ll{m^)+pl] 
1 

X 



Ll{m?) 



ht''^^\x,pl) = --ftT''^^\x,pl). 



(Bll) 



(B12) 



Dipolar form factor. The final results, already given in Eq. ((99|) and (jlOip . are 

f^q{a)( ^2^_9l M eg (1 - x)^ x{m + xM) 



hi^^^\x,pl) = --ftr''^^\x,pl) 



(B13) 



Exponential form factor 

2 Im J1 = -8P+X (m + xM) ga I^^ 
2P+ a; (m + xM) 



-.9a 



7rp^ 



2/^2> 



r 0, 



^a(»7^ 



2 Im Jl" = 8P+ (to + xAf ) ga IT^ 
2P+ (to + xM) 



9a 



■npi 



r 



{l-x)kl 
Llim') 



r 



'(l-x)A2 



r 



Llim^)+pl 
{l~x)Al 

il~x)Al 



(B14) 



The final result is, then. 



JlT [^iPt) ~ 



gl Mel x{m + xM) -[pl+Ll{^r^)]/[(i-x) hi] 



4 (2nfpl[Ll{m^)+pl, 



r(o,/"(-lVrro,"-(-^)+^- 



(l-x)A2; V ' (i-^)A2 



f-L<j(a),'_ ^2- 



/.^^^'^^(x,p^)^_±/^-L^«('^)(^,p2 



(B15) 



We have 



3. Axial-vector diquark including also longitudinal polarization 



.±,(a),.,9a{p^) 1 Mel 2ImJ« 

J IT \ ' "tJ ~ 



4 (27r)3 4(1 - x)P+ p2 _ rrfi 

,±,{a), OS ga(p') 1 Me^ 2ImJ{° 

1 ^ '^^^ 4 (2^)3 4(1 - x)F+ p2 _ ^2 



(B16) 



where the J° and J{" integrals are defined as in Eqs. (jB9p and (jBlOp . respectively, but now the second line in Eq. ([T 
is employed for the d^i,{p — I — P) and dua{P — p) Lorentz structures. 
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Point-like coupling (to avoid divergences we assume that the p^ integration is extended up to a finite cutoff A^) 
2 Im J1 =gaj -^ -^-^ {2m) SiD^) {-2ni) 6{D,) 



(27r) 
X Ix{x,pI) (1 + 

'It ■ P 



P 

W{x,pI) 



2MIJ -' y 2 

Yi{x,pI) +Y2{x,pl){lT-Pr)+Y^{x,pl)ll 



Ot ■ \i-T ^Pt\ } 1 



2ImJ[''^ga 



dH 1 



27r)4 Di D3 
X{x,pI) (1 + 

It ■ p 



{2m)5{D2){-'2m)5{Di) 



ll \ , ( _ 1 , 



Pt 



Yi{x,pi.) - Y2ix,p^) {It ■ Pt) - Y:i{x,p^^)^ 



+ 2Ml [(to + xM) (1 - x) - xm] + Aff [m {Ka - 1) - M(l + k^) cc] Z^ 



(BIT) 



where the integrals X,W,Yi, i = 1 — 3 are listed in App. O Unfortunately, most of the above combinations are 
divergent under the Mt integration. This is a typical pathology when choosing the point-like form factor for 
the nucleon-quark-diquark vertex, without any ad-hoc cut-off. 

Dipolar form factor 



2 Im Jf = gal -^ -^;i5^ (^TTz) 5{D2) {-2m) S{D,) 

Wix,pl) (^-^-^] St ■ [It - 2Pt] + 



pI 



Yi{x,pI) + Y2{x,pI) ^^^ + Y3{x,pI) ll 



Pt 



2ImJ('^ = g, 



dH 



1 



4 ^^^^^^2(27r*)5(Z?2)(-27r*)<5(I?4) 
xi -X{x,pI) (i + 3_]It-{pt-\It 






L T ' Pt 

p\ 



Y,{x,pl)-Y2ix,pl)lT-PT-Y3ix,pl)ll 



2Ml [(m + xM) (1 - x) - mx] + Ml [to {Ka - 1) - M{1 + «„) a;] l] 



(B18) 



where the integrals X,W,Yi, i ~ 1 — 3 are listed in App. [Cl Unfortunately, most of the above combinations are 
divergent under the dir integration, unless a dipolar form factor is considered with a higher degree, for example 
proportional to [p^ -KL^(A^)]~^ in Eq. P^ . This would introduce a 1/(1)3)'^ term inside Eq. (jBlSp . instead of 
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1/(1)3)^, and grant the convergence of the various integrals. With this very choice, we obtain 



2lmJ^=ga^-^-^lx{x,pl) 



2P+ 



2 \ -7-' dip 



2 7-/ dip 



2 \ -7-' dip 



( pI 1 

V2M2 2 






4M2 ■ 



7-/ rf2p 



+ 



Y^ix,p',)irP + Y,ix,p',)r,^'P + Y,ix,p'^)r, 



2 \ q-' dip 



j-i dip 



w{x,pi) x;"p - 2p, • 5,x; 



/ dip 



2lu,J['^ = g^\-^\x{x,pl) 



2P+ 



J2 \ -7-' dip 



yi(x,p^)X;"^ - Y,ix,pi)l',''^ - Y,{x,pi)X'^ 



-pIi'/'p 



2 \ -7-' dip 



1 



Pt 



2 2M„^ 



7-/ f/ip 



7-/ (iip 



4m: 



2-^5 



2 \ -7-' dip 



j-f dip 



2M^ [(m + xM) {l-x)~mx] X;"^ + M^ [m {na - 1) - M(l + k^) x] T 



f dip 



(B19) 



where the integrals X^' '^, i = 1 — 7, are listed in App. [C) 
The final result is. then, 



J IT \'^^Pt) 



9l Mel 



[l-x)^ 



32(2^)3 (p+)2[i2(A2)+p2]3 



X{x,pI) 



p2x;*p + 



/ Pr _ ]\ j'dtp ]_ 

{2M^ 2) 2 4M2 



yi(x,p2)i;*P + Y2{x,pl)l'f^ + y3(a;,p^)X^*^ + W{x,pI) (x;** - 2p, . S^I', 



T-i dip 



dip 



ht'^^\x,pl)^~^^^' 



{1-xY 



32(2^)3 (p+)2[^2(A2)+p2]3 



X(a;,p^) 



2 -7-/ dip 
Pt-^1 



1 p? 



2 2Af 



/f2 



7-/ dip 



1 



7-/ dip 



4M2 5 



+ Yi{x,pI)I['''p - Y2{x,pI)I'/'p - Y3ix,pl)l'/'P + 2M^ [(m + xM) (l - x) - mx]l['''P 
+ Ml [m {Ka - 1) - M(l + Ka) x] I'/'P \ . (B20) 



Exponential form factor 

2 Im J1 = ga / 7—4 ^^^ (2^i) <5(i?2) (~2^*) (5(i?4) 



(2^)4 
X \x{x,pI) (1 






jTT/ 2 \ / 't ' '-'t 



*^T ' I't '-'PtI 



I'jp • J)rj 



Pi 



Yi{x,pI) + Ij. ■ PtY2{x,pI) + IIY3{x,pI 



2P^ 



x{x,pi 



2 q-exp 1^ / Pt \ q-exp ^ q-exp 



ri(x,p^)i-p + y2(x,p^)xrp + r3(x,p2)i-p. 



2IniJ("=g, 



W{x,pI) {lT''-2p^-S^Il 

dH g[(p-/)'-m']/A: 



{2^Y 



D1D3 



{2TTi) 5{D2) i-2T:i) S{Di) 
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x|-^(-,P^)(l+2i^)^.-(p.-^^. 



t,rj-i ' ^-Y 



Yx{x,p\) ~ It ■ PtY2{x,pI) - llY:i{x,ply 



+ 2Ml [(m + xM) (1 - x) - mx] + Mj [m{K.a - 1) - M(l + k^) x] ll 



2P+ 



X{x,pI) 



2 q-exp \_ ( Pt \ 'j-exp ^^ ^ ^-rexp 



Y,{x,pI)1T^ -Y2{x,pl)ir^ -Y,{x,pI)IT''+ 

2Ml [{m + xM) (1 - x) - mx] ll'^^ + M^ [m{Ka - 1) - Af (1 + «„) i'] ^a'"'' 



(B21) 



where the integrals X, VF, Fi, i = 1 — 3, and X^^" 
The final result is, then, 



1 — 7, are listed in App. [Cl 



pJ_q(o) 



5a 



AV^"^(a:,p^,) = -^ 



2 Me^ e~[PT+ia(™')]/[(i-^)A'] 



32(2^)3 (p+)2[i2(^2)+p2] 



X(a:,p2) 



,2 Texp 1^ / Pt 



Pt-'-i 



V2M2 



- xr^ 



1 



4M2 5 



+ yi(x,p^)i-p + r2(x,p^)ir'' + Y,{x,pl)irp + w{x,pI) {IT" - 2p. • s^xf') 






2 q-exp 






^ Pt \ q-exp 1^ ^ q-exp 



32 (27r)3 (P+)2[L2(TO2)+p2] 

+ yi(x,p2 )i-p _ r2(x,p2)i-p _ r3(a;,p2)i-p ^ 2^^2 j(^ _^ ^^^) (1 - x) - mxixr" 



Ml [m{Ka - 1) - Af (1 + Ka) x] 1^ 



(B22) 



4. Axial-vector diquark including also time-like polarization 



We have 



ht'^''\x,pl) 



ga{p^) 1 Mel 2lmJf 

4 (27r)3 4(1 - x)P+ p2 _ ^2 

5a (p^) 1 Mel 2ImJ{° 



(27r)3 4(1 - x)P+ p2 „ ^2 



(B23) 



where the Jf and J(° integrals are defined as in Eqs. (|B9p and (jB10[) . respectively, but now the last line in Eq. PH]) ') 
is employed for the d^i/(p — I — P) and d^a{P — p) Lorentz structures. 

• Point- like coupling (to avoid divergences we assume that the p^. integration is extended up to a finite cutoff A^) 

2ImJf = -4:P+x [m(2K„ + l) + Af (2Kaa; + l)] galf' 

P+ X [m {2Ka + 1) + Af {2Ka X + I)] , f Ll{m^) + pi 
= -.9a n log 



Trpj 



■p.i. 



2ImJ('' =4P+ [m [(2Ka - 1) a; + 2] + xAf [(Ka - l)2x + 3]] galf 

P+ [m[(2Ka-l)a; + 2]+a;Af [(Ka - 1)2j; + 3]] f Ll{m'^) + pi 



np^ 



Llim^) 



(B24) 



where I^' ' is the same integral as in Eq. (|B4p but with the substitution Lg{m?) ^-> Lairn^)- 
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Using again Eq. ([8]), the final result is 

J IT K^iPt) — ~ /r,_\4 ZTTTTTZTTTTTrTi ^"S 



(2^)^ 



p2 [L2(m2)+p2] 



Ll{m^) 



,-Lg(a) 



(a;,p2) 



?2 Afe^ TO [(2k„ - 1) a; + 2] + xAf [(k„ - 1) 2a; + 3] , /L2(m2)+p; 



(27r) 



vllLlim'^) + pI 



log 



L2(to2) 



(B25) 



• Dipolar form factor 



2 Im J° = 4P+a; (1 - x) [to (2k„ + 1) + M (2Ka a; + 1)] 5a if ^ 

__ P+ x{l-x) [to (2Ka + 1) + Af {2Ka X + 1)] 

'•'^" ^i2(A2)[i2(A2)+p2] 

2ImJ{° = -4P+(l-a;) [m [(2k„ - 1) a; + 2] + xM [(«„ - l)2a; + 3]] fla^f^ 
P+ (1 - x) [to [(2Ka - 1) .X + 2] + xM [(«„ - 1) 2a; + 3]] 



-ffa' 



^i2(A2)[i2(A2)+p2 



(B26) 



The final result is 



f^ 9(a) ^ 2 N 3a A^ e? a: (1 - a;)2 ["^ (2Ka + 1) + M (2k„ a; + 1)] 



I IT 



8 (27r)4 P2(A2)[i2(A2)+p2]3 

.9^ Af el (1 -- x)2 [to [{2Ka - 1) .T + 2] + xM [(«-,„ - 1) 2a; + 3]] 



8 {2ttY 



L2(A2)[L2(A2)+p2]3 



(B27) 



We find a discrepancy between these results and those of Eqs. (18) and (24) in Ref. [46], probably due to errors 
in that calculation. 



• Exponential form factor 

2 Im J{' = -4P+.X [to {2Ka + 1) + M {2Ka x + 1)] ga IT^ 



= -9a ■ 



P+ X [to {2Ka + 1) + M {2Ka X + 1)] 



t^Pt 



r 



Ll{m^) 
(l-x)A2 



-r 



iaK)+P? 



(l-a;)A2 



2Im J{°- = 4P+ [to [{2Ka - 1) a; + 2] + a;Af [(k^ - 1) 2a; + 3]] ga^T 



9a 



P+ [to [{2Ka - 1) .T + 2] + a;Af [(k^ - 1) 2a; + 3]] 



■Kp- 



r 0, 



(l-a;)A2 



-r 0, 



Ll{m')+pl 
(l-x)A2 



(B28) 



The final result is, then. 



f±q(a), 2^ ^ 9l Mel x[m{2Ka + l) + M{2KaX+l)] ^^^2^ + ^2j^2^y^^^^^^f^2j 

hT [^,Pt) g (2^)4 p2[i2(^2)+p2] 



2^^2^ 



r 



Li{m 



r 0, 



(l-a;)A2 
l]+x^ 

Llim') 



r 



,±q(a)f 2^_ 9^ M el ui [{2Ka - 1) a; + 2] + a;A// [{Kg - 1) 2a; + 3] -[pi+L^(m-')]/[{i~x)A^] 

' ^'''^^' 8 (2^)4 p2[i2(^2)+p2] 



(l-a;)A2 



-r 



Llim^)+Pl 
{l-x)Al 

Llim^)+pl 
{l-x)Al 



(B29) 



APPENDIX C: USEFUL INTEGRALS 



In this appendix, we calculate the relevant integrals that repeatedly show up in the expressions of T-odd parton 
densities for all choices of nucleon-quark-diquark form factors, when vector diquark propagators are represented in 
the first and last forms of Eq. (fTO| . 
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We will systematically use the substitutions l'^ ^ Ij. — Pt and y = l'^ + L^(m^), or y = l'^ + L\-{K\) for the 
dipolar form factor, where Lx is defined in Eq. ([5]) for X = s^a scalar and axial-vector diquarks, respectively. We 
will also encounter the following angular integrals, where 6 is defined as the angle between l'^ and p^, and (/>, 0s, are 
the azimuthal angles of p^ and St with respect to the scattering plane: 



de- 



2ir 



de 



I^^IIPtI cos 6* + p^ 
VT^+pl + 2\l'^\\pT\ci 

[\V^\\pT\cose+pl] 



7r(l-sgn(|Z;|-|p^|)) , 



I' 



Pt 



2|i;,||pr|cos( 



- H;2 + 3p2 



f 

Jo 



'" d0 l^^ll-^^ l ^o^[^ + ("^ - ^^^^ + \Pt\\^t\ cos( 



65) 



/2 



I' 



Pt 



2|i;,||p^|cos( 



Pt 



\St\ , 

n cos( 



IPt 



(/.s)(sgn(|p,|-|Z; 



^" (IZ^IIg^l cos[g + (0 - 0s)] + |p^||5^| cos(0 - fe))- 
Z;2+p2+2|Z^||p^|cos0 



2p2 



p^-Z^^ + |p^-Z^2|) cos 2 



2p^ 



(CI) 



Point-like coupling 



j-p.i. 



dl'^ {1't+Pt)-P7 



(27r)2 
1 



Pt 



{Vt+PtYW + L\{iv?)] 



{2^)'pI Jo 



'^''^"^^'z^2 + L2,(m2)y^ 



i'2^)PlJo 



a|ty ||ty| 



1 



d0 



|Z;||p^|cos6>-fp2 
Zi,2+p2 _^2|Zi,||p^|cos( 



Lx{m')+PT 



dy 



47rp2 Jl2^(^2) 



log 



47rp2 



L2^(m2)+p2 



(C2) 



Dipolar form factor 



-^1 



dZ; (Z; + P:j 



Pt 



(27r) 



Pt 



[l'T+PT?[l'T^+Ll{K\)f 



LliA-D+p^T 



(1-x) /"^-^'^''xl+P-T dy 

47rp2 /,o ,.0 , 1/2 



(1-x) 



4^i2^(A2,)[L2^(A2,)+p2] 



(C3) 



Exponential form factor 



^exp 



dl'T {1't+Pt)-Pt e-['^'+^i('"')l/[(i-)Ai] 



(2^)2 
J_ /-Ipt 



p2 



(27r)p2 



d\Q\l'^ 



{I't + Pt? [I't^ + Llim^) 

p^ll'T^+Ll{m^)]/lil-x)Al] 



T Jo 



Z;2 + L2,(m2) 



(47r) p2 



T JLj,(m^) 



y 



47rp2 



r 



£2^(m2) 
(l-^-)A| 



-r 



L2^(m2)+p2 
(l-x)A' 



X 



(C4) 



Next, we list the coefficients and calculate the relevant integrals that are needed to construct T-odd parton densities 
for all choices of nucleon-quark-diquark form factors, when vector diquarks are represented in the second form in 

Eq. uni). 
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X^- 



w 



Yi 



8(m + M) {P+f (Ka - 1) (1 - x) 



Ae 



8(P 



+ ^2 



M2 



Yo = - 



Ka hm^ - 3xm^ + Mx^rn^ + 2M7r? - iMxru^ - M^xm + 2M^xm + M^x^ - M^x 

- [pI + Ll{m^)] (M + m)) + (1 + na)Mlx [M(l +x) + m] + [pi + Ll{m^)] m (k„ - 1) {x - 1) 

+ m (Ka - 1) x^ {Ml + m^) + m^ {2x - 1) + M^m (1 - 2a; + (1 + Ka)x'^) 
8m{P+f{Ka-l)il-x) 



M2 



Y. 



4{P 



+ ^2 



Mf 



Kq f - m^ + a;m^ - Mx^rr? - 2Mm^ + iMxm^ - M^m - 2M'^x^m + iM'^xm ~ 2Mlxm 
- M^x^ + [pI + Ll{m^)] (m + M) + M^x) + Mlm(Ka - 1) - MMl{l + Ka)x , 



'T-/ dip 

^1 ^ 



dZ;, (Zi,+PT)-PT 



(2^)2 

-d\i',\K\ 



Pt 



{l'^+p^)^l'^^+Ll{Al)] 



2^13 



1 



(2^)2p2[Z;2+i2(A2)]3 



27r 



^6* 



P-ld|Z^||Z^| 



1 



1 



\l 'J\pr\cOSe + pl 

-pl + 2\l'^\\pj.\cost 
dy 



(27r)p2 [i/,2 + L2(A2)]3 A^pUl^kd y 



Ll(f^l)+Pl 



SttpI 



2Ll{Al)+pl 



87rL4(A2)[i2(A2)+p2 



(C5) 



^2 ^ 



dl' 



{I'r+PT? 



(2^)2(Z;+p,)2[Z;,2 

{2^Y %^ + Ll{hl)YJ, 






de ^ — 

An 



dy 



&^Lt{^l) 



(C6) 



'T-/ dip 



dX [X+PT)-Pr 



1 



(27r)2 
"= rf|Z^||Z^ 



Pt 



(Z^+P,)2[Z;,2+L2(A2)p 



(27r)2p2[i^2+i2(A2)p 



de 



\V^\\Pt\ COS 9 + pI 



1'^^+pI + 2\1'^\\Pt\cos0 



1 



Anp'^ 



\pt\ 



d\l'J\l', 



[I? + LI 



2(A2)]3^^P- i.) + P./ ^^, 



d|z^||z;i 



.^ + i^(A2)]3 
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87rp2 



r-ia(Aa)+P? 



I Ll(Al) 



dy 



2pl + LliAl)-y 



2p^ + SLl{Al)p^ 
S7rpl\\2Ll{Al)[LliAl)+pl]^ 

Ll{Al) + 2pl 
167rL4(A2)[L2(A2)+p2]' 



P. 



P 



1 



nmAD+p 



.212 



J-/ dip 



dl'^ {1't+Pt)-Pi 



{I't+Pt 



(2^)2 

-d\v,\K\ 



Pj 



{l'^+p^YW + Ll{Al)Y 



{2i:YpIW + LI 
d\l'^\\l'^\ 1 



1 /■ / 

^2(^2)13 7 d^(|i;i|PT|cos0+p^ 



(2^) [Zi,2 + L2(A2)]= 



47r 



dy 

L2(A2) y^ 



1 



8^i4a(AD 



'T-/ dip 



7-/ dip 



dl' 



{K+PtY 



{2i:f{V^+p^Y[V^^ + Ll{Al)f 

■ 2-K 



T T 



1 



(2^)2 %^^Ll{Al)f}^ 

d\V^\\V^\ {I'^^+pD 
(27r) [Zi,2 + L2(A2)]3 



8^if(A2) ' 



^6* U 



1 

47r 



-p?, + 2|Z;||pr|cos0 
y-i2(A2)+p2 



dy 



L2(A2 



W dip 



< (Z^+p^)-5x 



1 



(27r)2 S2 



(Z^+p,)2[i;2+^2(A2)]3 



1 



i2nrpl%^+Ll{AlWJo 
1 |5.| ,, _ /-^^^(^"Hp- dy 

2L2(A2)+p2 



'" ^0 I^tII^t I COs[0 + (0 - (/-s)] + \Pt\\S^\ COs(0 - (/)5) 



l'^^+pl + 2\impj.\cose 



1 



SttS'? \p^ 



S-kSI 



[Pt ■ St 



L4(A2)[L2(A2)+p2]2' 



LtiAl) [i2(A2)+p2 



COS 



7-/ dip 



dl'^ [{I't + Pt) ■ St] 



1 



(i;+p,)2[i;2+^2(A2)]3 



1 



de 



\1't\\St\cos[9 



(2^)2 ^ 

(27r)2S2 [Z;2 + i2(A2)]3 7^ - 1',^+pI + 2\1',\\pt\cos9 

1 I /-i(A^)+p| dy 52 

L2(A2) 



|Pt||S't|cOs(0-(/)5) 



t2 




„ gTT^ [p^+(p2+i2(A2)_y)cos2(0-</.5)] + 7r52 
^ y Pt 



dy 



22 

T "T 



52 pi 



4L4(A2)[L2(A2)+p2] 
^4[L2(A2)+p2]2 



i2(Al)+PT ^^^ 

2 ^ (cos 2(,^ -^s) + 1)pI + (cos 2((/) - (/-s) + 2) i^ (A2)) 



1 



167rL4(A2)[L2(A2)+p2 

1 

167rL4(A2)[L2(A2)+p2 



L 






LiiK 



iJrj-, 
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(C7) 



q-exp 



'f-GXp 
^1 



f dl'^ e-l'r+-f^a(™')]/[(l-==)Aa] 

1 /-oo _ e"^/[(i"=^)^'l 1 



47r JLlim^) 



V 



i^nO' 



47r 






1 ' ' ' 






47rp| 



!y IPtI COS0 



P. 






Tit 2\ I 2 



StTP?. [iL2(m2) 






ia(™^)+P|- 



dy 



g-y/[(l-x)A^] 



1 J'rr2^ 2a , o„2iT^/^n -^a("^^) \ r.2. 2a , „2^j.fn Lljm^) + pI \ 



q-exp 



< e-[J^^+L^„(m^)]/[(i-.)A^] [{V^+p^).p^ 



^exp 



(27r)2 (i^+p^)2[Z/,2+i2(m2)] 
(27r)2p2 [i/,2+^2(^2)] 



Pj 



[I't + Pt 



27r 



de (|z;i|p^|cos0 + p? 





1 

47r 



2^ [Z;2+^2(^2)] 

g-y/[(l-x)A2] -^ 



dy- 



L2(m2) 



y 



£2/ 2) 

r — — — — 

Att V (1 - a;)A2 



q-exp 
-^2 



dZ; e-liT+Ll(m^)]/[(i-^)f^l] 



(2^ IzT+p^m^+^fR)! ^^" ^ ^"^ 

|Z;||ZL| e-['T+i^(™^)]/[(i-)A^] /•2- 



io (2^)^ [Z^2+i2(^2)] y^ I: 

in J r 2(^2^ " y" 



r/2 

T 



-P. 



-2|Z^||pr|cosf 



'i2(m^) 



'T-r-aV"'')]/[(l-:r)AS] 

[Z^2+i2(^2)] ^T +PtJ 



41 



q-exp 
-^6 






(2^)252 [l'^^ + Ll{m^)] Jo l'^^+pl + 2\imp^\cose 



Ll(m^)+pl p-y/[{l-x)Kl] 



Att sIpI JLKrn^) y 

Pt-S^ /rfo^ ^^K) \ j.(QLl{m^)+Pl 



X 






■e2;p 

d|Z;||i^| e-['i-'+^^('"')l/[(i-)A^l [^^ ^^ \l',\\Sr\cos{e + (/. - 0s) + Ip^IIS-^I cos(0 - ci>s) 



(2^)252 [z;2 + ^2(^2)] j^ v^^ + pi + 2\i'^\\p^\cose 

1 f /•ia("l^)+PT e-2'/[(l-^)^'l ^2 

(27r)2S'2 [7l2(„2) y Pt 

/•°° g-y/[(l-x)A^] ^ 

+ ttS^ I dy > 

JLl(m^)+pl. y J 

87rp2 1 V J 

+ r (0, ^ii^) [[p^ + L2 (^2)] cos 2(0 - 05) + P^] 

'^(°' i!'"l)lf ) [P^+ia(m^)] COS 2(0-05)1. (C8) 
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